“Calhoun 


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1965 


Numerical examples in the investigation of a 
particular matrix in eigenvector theory. 


Yost, Albert N. 


Monterey, California ; U.S. Naval Postgraduate School 
http://ndl.handle.net/10945/12740 


This publication is a work of the U.S. Government as defined in Title 17, United 
States Code, Section 101. Copyright protection is not available for this work in the 
United States. 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 
| (8 D U DLEY research materials and institutional publications created by the NPS community. 
«ist sia Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


NY KNOX appointed — and published -- scholarly author. 

ia) LIBRARY Dudley Knox Library / Naval Postgraduate School 

411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 


NPS ARCHIVE 


1965 
YOST, A. 





NUMERICAL EXAMPLES IN TRE INVESTIGATION OF 
A PARTICULAR, MATRIX IN FIPENVECTOR THEORY 


“ALBERT N, YOST 


DUDLEY KNOX LIBRARY 
NAVAL POSTGRADUATE SCHOOL 
MONTEREY CA 93943-5101 

















NUMERICAL EXAMPLES IN THE INVESTIGATION OF A 
PARTICULAR MATRIX IN EIGENVECTOR THEORY 


ewe oR OK OR 


Albert N. Yost 





NUMERICAL EXAMPLES IN THE INVESTIGATION OF A 


PARTICULAR MATRIX IN EIGENVECTOR THEORY 


by 


Albert N. Yost 


Lieutenant, United States Navy 


Submitted in partial fulfillment of 
the requirements for the degree of 


MASTER OF SCIENCE 
IN 
AERONAUTICAL ENGINEERING 


United States Naval Postgraduate School 


Monterey, California 


L Oet6n 





NUMERICAL EXAMPLES IN THE INVESTIGATION OF A 
PARTICULAR MATRIX IN EIGENVECTOR THEORY 
by 


Albert N. Yost 


This work is accepted as fulfilling 
the thesis requirements for the degree of 
MASTER OF SCIENCE 
IN 
AERONAUTICAL ENGINEERING 
from the 


United States Naval Postgraduate School 





ABSTRACT 

For every orthonormal matrix X there exists a skew~symmetric 
matrix N such that X = (t-n)(14+n) |, provided that (r+n5 + exists. 
A matrix K, similar to N, can be defined for biorthonormal matrices 
U and V such that U = (r-k) (1+K) +, provided that (1+K) "+ exists. 
Numerical methods are presented for examination of the properties 
of K. The particular property anticipated for K, that it exhibit 
n(n-1) basic parameters inherent in biorthonormal matrices, is not 


apparent in the numerical examples derived. 
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1. INTRODUCTION 


If A is a real symmetric matrix of order h, it possesses n real ortho- 
gonal eigenvectors x, which may be norm&lized, as described in Appendix 1. 


The square modal matrix X' of these eigenvectors, 


' — 
X (x) Ky sors x) (1) 
is orthonormal, that is: 
X'x' = 1 =x! X (24 


For the purpose of this analysis the eigenvalues of A are assumed to be 
discrete and the eigenvectors unique. 

The Ac elements of such an orthonormal matrix are subject to the 
conditions of equation (2), of which only n(n+tl)/2 are independent, leaving 
n(n-1)/2 independent parameters by which the orthonormal matrix can be ex- 
pressed. Consequently, as has been shown by Heading [3] for example, X' 
can be written in general as 

x' = (1-N) (14n) 7+, (3) 
provided that (1+n)} exists, where N is skew-symmetric and contains the 
required n(n-1)/2 parameters. Hence any real skew-symmetric matrix de- 
fines an orthonormal matrix, provided that (I+N) is non-singular. Equa- 
tion (3) may be solved for N(X"): 

N= (14x') Tex") (4) 

The proof of equation (3) hinges upon the fact that I+N and I-N 


commute, i.e. 


~ =] =] 
re (I-N) ~(I+N) (L-N) (I+N) 


(r-n)*(1-N) (14N) (14) 72 

=e 
Since (I+N)(I-N) = Tne = (I°-N)(It+N); furthermore N must be skew-symmetric 
in order that I-N = I+N as used above. 


The discussion of this property of orthonormality in the context of 





the eigenvalue problem is of particular significance in the reversed prob- 
lem suggested by Bell [4]. The solution of his problem for A symmetric is 
made possible by the existence of N, as in equation (3), containing just 
the right count of parameters required for the solution. As pointed out 
in [4], the n(ntl1)/2 independent elements of a symmetric A should yield 
the same total number of independent parameters in the eigenvectors and 
eigenvalues of A. This exact count does appear in the n eigenvalues plus 
the n(n-1)/2 elements of N. 

This description of the property N(X'), or N(A) for A symmetric as 
the property is restated in the preceding sentence, is background for the 
work to be reported herein. In [4] Bell has shown that an attempt to ex- 
tend his analysis to a class of unsymmetric matrix A, possessing a real 
eigenvector matrix U', would involve the property K(U') analogous to N(X'), 
where the matrix K should express n(n-l) independent parameters inherent 
in U', or twice as many as are inherent in X'. The development of this 
concept is presented in the next section, where the multiplicity possible 
for K is pointed out, as is the need for a definition, if possible, of 
that canonical form, K.> which displays the required property. 

Numerical methods are set up in the present investigation for de- 
riving K(A), with A unsymmetric, and for exploring its properties in 
association with those of the modal matrix U', and of the modal matrix 
V' which is adjoint to U'. These methods are programmed for arbitrary 


matrices of order up to n = 25. 


2. PROBLEM STATEMENT 
If A is real but unsymmetric, then its eigenvectors u, are not 


orthogonal; but the adjoint set of vectors Va? belonging to A, satisfies 





a biorthogonality with the u;, as described in Appendix 2. In fact the 


Square matrices of these adjoint vectors, 


1s 
U (uy) Up ees u_) 
| (5) 
V =(v, Vo so Vn) 
satisfy the conditions of biorthonormality: 
NA 
Viul = 1 =u'v',” (6) 
if the us and v, are mutually normalized. 


It should be noted here that, as before, the eigenvalues of A are 
assumed to be discrete and the eigenvectors unique; A is also restricted 
to that class for which the eigenvectors are real. 

The 2n* elements of such real adjotne matrices in (5) are subject 
to the x2 conditions of equation (6), all of which are independent, indi- 
cating that both U' and V' can be expressed with no more than n* basic 


parameters. As shown in [4], U' and V' can in fact be written as 


U' = (1-K) (14K) > 


a (7) 
v' = (I+K) (1-K) 


provided (14x)! exists, where K is an arbitrary real matrix. Hence , by 
equations (7) any real matrix defines a pair of real adjoint or biortho- 
normal matrices U' and V' under the provisions stated. Equations (7) 


may be solved for K: 


K = (rtu') b(r-u') 


zs (8) 
& = (v'41) lqv'-r) 


It is readily verified that equations (7) and (8), for A symmetric, 


are identical with equations (3) and (4), since for A symmetric: 


Dp = art = pa 
oS co) 

K = N #2 -K 
The proof of equations (7) follows from the commutability of (I+K) and 
(I-K), as for equation (3). The consistency of equations (8) is demon- 
strated for the general case in [4]. 

The diagonal elements of I in equation (6) are the normalized inner 

product, or dot product, of the gett columns, or vectors, of U' and V'. 


Unlike the case for the vectors x, in equation (2), this is not a unique 


i 
specification for the dot product of us into itself, or of Vs into itself, 
i.e., for the "magnitudes" of uy and V.- For each such pair of adjoint 
vectors u, and ae (belonging to the same eigenvalue of A) there is a free 
condition to be chosen arbitrarily. There are n such pairs and hence n 
such arbitrary conditions. This suggests that the property K(U') can in- 
volve not as many as na basic parameters, but rather as few as n(n-1). 
An additional argument given in [4], on the basis of the eigenvalue problem 
itself, is similar to the one given in the Introduction for symmetric A, 
namely: an unsymmetric matrix A has ae independent elements and should 
yield just as many total basic parameters in its eigenvalues, which are n 
in number, and in its eigenvectors, which must therefore have only oe 

As shown in Appendix 2, there are an infinite number of pairs of 
biorthonormal matrices, satisfying equation (6), belonging to a given 
matrix A. This multiplicity corresponds to the fact that the eigenvectors 
can be multiplied by any scalar and still satisfy the eigenproblem state- 
ment, coupled with the free condition on the dot product mentioned in the 


preceding paragraph. In view of equations (8), there must be a correspond- 


ing multiplicity in K(U') for any A, or in other words in K(A). 





If K can be found to involve as few as n(n-l) independent parameters, 
the skew-symmetric form which K takes with A symmetric suggests that in 
the more general case these n(n-1) parameters would be off-diagonal elements, 
and that the main diagonal would remain empty. This form of K, if it should 
occur only for one particular pair of U' and V' out of the infinite multi- 
plicity possible, might be designated its canonical form, Ko» corresponding 
to the canonical pair, Ul and Vie The conditions on U' and v' for the 
canonical form, if it exists, have not yielded easily to analysis for n 24, 
but might become apparent from a study of numerical examples. 

The problem which was stated for this investigation may now be outlined 
as follows: 

1. Set up the computer solution for K from a given unsymmetric A, 
using equations (8). 

This has involved the use of known techniques for diagonalizing 

A; Laguerre's method [51 was adapted. It was also necessary to de- 

velop subroutines. for numerical solutions of the adjoint vectors U' 

v', when the eigenvalues by Laguerre's method were known. Care must 

be taken that the given A possess real eigenvectors. The solutions 

for U' and v' are verified for biorthonormality in equation (6). 

2. Set up the computer solution for U' and v' from arbitrary K, 
using equations (7). 

This is provided as a final subroutine to Step 1, for an over- 

all check of the solution; it is also provided as a separate program 

for the investigation of U' and ve from arbitrary K.| The latter pro- 

gram also includes solutions for V'v', U'u', v'¥',u"U! and V'u' (see 

Appendix 6) to provide for possible interpretations of the numerical 

relationships among the adjoint vectors themselves. 

The objectives of Step 1 were first to provide an overall check of 
the methods developed in [4] for A unsymmetric, and second to exhibit 
numerically some examples of the property K(A), with particular regard for 


evidence of its basic parameters. The objectives of Step 2 were to explore 


whether evidence of the conditions on these parameters might be exhibited 





readily in the numerical characteristics of U' and V', for example in or 
among the ''magnitudes" of the eigenvectors, as displayed in the elements 
of the diagonals of U'u' and V'v'. The computer programs for these in- 
vestigations are written in Fortran 63 and can accommodate matrix orders 
of 25 or less. 

The matrix properties which are being sought out in this investiga= 
tion are known to have hypergeometric interpretations for n greater than 
three [4]. For n less than three they represent the familiar interpreta- 
tions of Euclidean three-space. Attempts to extrapolate the latter into 
higher order hyperspace can be misleading. On the other hand, it was con= 
sidered desirable to keep n small at the outset in order to limit numbers 
of numerical data to be examined. For this reason the numerical examples 
reported herein are for third and fourth order matrices only. 

A third computer program was formulated in order to further investig- 
ate the characteristics of the matrix K. In this case the K matrix, how=- 
ever, was derived from an orthonormal matrix after the latter had been 
premultiplied by a diagonal matrix. The goal here was to observe the ef- 
fects of premultiplication in the changes on the matrix K. Appendix 3 


describes in detail the background for this program. 


3. COMPUTER PROGRAMS 
Two Fortran 63 computer programs were formulated in order to satisfy 
the requirements set forth in the Problem Statement above, Steps 1 and 2. 
All programs were arranged to operate with matrices of order n < 25. 
Appendix 5 describes the computer program named REVEIG, designed for 
Step 1, which basically is a general eigenvalue program for any real square 


matrix A where, as described in Appendix 2, 





Au' = u'D 


Av' = v'D 
and Yu! eer 


The derived matrices U' and V' are then applied in equations (8) to derive 
K and ‘kK, from U' and v’ respectively. 

The available EIG3 subroutine used to derive the eigenvalues for REVEIG 
can be used to find real and imaginary roots of a matrix and store them for 
further use [5]. A basic assumption for REVEIG was that all matrices in- 
vestigated would be required to have real roots only. Subroutine EIG3, 
however, did not derive the eigenvectors along with the eigenvalues. A 
program for the eigenvectors was derived as an additional step, as described 
in Appendix 5. It should be noted here that due to the arbitrariness of the 
solution used for the eigenvectors, a particular restricted class of multip- 
licity is built into the program, to the exclusion of other classes which 
exist, and which might bear further investigation. Some typical test re- 
sults from REVEIG are contained in Appendix 8. 

Appendix 6 describes the computer program called REVEIG1, designed 
for Step 2, which was based on equations (7) in the Introduction. Inputs 
were made to REVEIG] from the K and K matrices derived in REVEIG, in 
order to provide a check of REVEIG1 by recomputing the U' and v' matrices 
from which K was derived. Thereafter, other arbitrary matrices could be 
used as K inputs to REVEIG]. As a part of this program the products V'v', 
v'y' : Tu! : u't! , and V'U' are also derived. Some typical results are 
shown in Appendix 8. 

Appendix 7 describes the computer program called REVEIG2. This 
program was formulated late in the investigation in order to further explore 


the effects of matrix algebra on an orthonormal matrix. As described in 





Appendix 7, the orthonormal matrik U' derived in REVEIG was premultiplied 
by a diagonal matrix G of significance in a particular case. The result 
was used as an input to derive the matrix K using equation (8) in the 
Introduction. The ultimate goal was to investigate the effects of pre- 
multiplying an orthonormal matrix in this special case, and to determine 
what characteristics, if any, exist in the resulting K matrix. The signifi- 
cance of this case is discussed in Appendix 3. 

The three computer programs are designed for an input format of 
(4E20.10). The outputs are all printed out in a (7E17.10) format. In order 
that the computer results in Appendix 8 be more easily understood, the fol- 


lowing relationships are defined: 


Computer Theory 
UP U' 
VP pS 
KT 

UPT vw 
VPT v' 


Comment cards are included in the program decks for all the programs, that 
basically describe the individual steps involved. The program decks can be 
acquired from the Department of Aeronautics of the U. S. Naval Postgraduate 


School. 


4. RESULTS AND CONCLUSIONS 

A few of the several numerical examples run from the programs described 
are included in Appendix 8, showing typical results. It was a primary obt 
jective of this investigation to provide an overall numerical check of the 
methods developed in [4] for A unsymmetric, namely to derive a property 
K(A) as described in Section 2. In this objective the investigation is 
successful, and the consistency of equations (7) and of equations (8) is 


verified. 





A second objective was to exhibit typical examples of the matrix K 
with the view that typical patterns might be evident in its properties. 

In particular the so-called canonical form dependent upon n(n-1) basic 
parameters was of interest. This form was not observed, at least not in 
the form with empty main diagonal, and was not otherwise identified. In 
this connection it should be noted that a single restricted class of multi- 
plicity has been built into REVEIG, as discussed in Section 3. 

Inputs were made into REVEIG1 consisting of arbitrarily assumed 
"canonical" forms, K.> possessing zero main diagonal. It was anticipated 
that the final printout of this program, consisting of the several products 
listed in equation (2-8) of Appendix 2, might show evidence of associated 
canonical properties in U' and V'. These properties were not observed. 

Finally a calculation with REVEIG2 was made to explore the form of K 
which results from a special case of unsymmetric matrix, as described in 


Appendix 3. No significance was found in the resulting example. 





BIBLIOGRAPHY 


Crandall, S. H. Engineering Analysis. McGraw-Hill, 1956. 
Lanczos, C. Applied Analysis. Prentice-Hall, Inc. 1956. 


Heading, J. Matrix Theory for Physicists. Longmans, Green and Co., Ltd. 
1958. 


North American Aviation, Inc., Los Angeles Division. A Feasibility 

Study on the Generation or Prediction of Structural Influence Coefficients 
to Meet Simple Aeroelastic Criteria, by R. W. Bell. January, 1964. 
Report No. NA-63-1183. 


Beresford Parlet, Mathematics of Computation, American Mathematical 
Society, Laguerre's Method Applied to the Matrix Eigenvalue Problems, 
Vol. 18, No. 87, July 1964: 464-485. 


10 





NOMENCLATURE 


Most terms within this thesis are defined when used; however, the 


following list of the general terms used is supplied for ready reference. 


X* 


tt 


real square matrix 
elements of A 

B’ ap"? = a derived matrix in Appendix 3 

diagonal matrix of non-vanishing, real diagonal elements 
diagonal matrix of eigenvalues assumed to be distinct 
diagonal matrix 

diagonal matrix 

diagonal matrix 

identity matrix 

matrix property belonging to A unsymmetric 

hypothetical canonical form of K 

integer 

skew-symmetric matrix 

eigenvector (column matrix) 

modal square matrix of eigenvectors u, 

hypothetical canonical form of U 
eigenvector (column matrix) 

modal square matrix of eigenvectors vy 
hypothetical canonical form of V 

diagonal matrix 

eigenvector (column matrix) 

orthogonal modal matrix of eigenvectors x 


af 
pseudomodal matrix (Appendix 3) 


ll 





real, discrete eigenvalues of A 
transpose of ( ) 


normalized form of matrix ( ) 


TZ 
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APPENDIX 1 


SIMPLE EIGENVALUE PROBLEM WITH A SYMMETRIC 


The simplest possible eigenvalue problem is defined [1,21 by 


AX} - AqX} (lel) 
where (see Appendix 4 for definition of terms): 


A =n x n Symmetric matrix 


eigenvalue for which a solution x; exists; there are n 
solutions assumed to be discrete 


hi 


et column matrix (eigenvector) 


The n solutions may be combined into a single statement 


AX ‘= xD (1-2) 
where 
X = modal matrix of eigenvectors 
D = spectral matrix (diagonal) of eigenvalues 


It is to be noted here that post-multiplication of the modal matrix 
by a diagonal matrix does not alter the eigenproblem (1-2). 
~ 
After premultiplying equation (l-2) by X, and also after trans= 
posing equation (1-2) and postmultiplying by xX: 
ad A/ 
XAX = XXD 


DXX 


Putting XX = W and subtracting: WD - DW = 0, i.e. Wij = 0, r5, 


since the )j are assumed discrete. Therefore 


XX = W (diagonal) 


a 
i.e., X is orthogonal. Letting xX! = XW 7, i.e., after normalization: 


ay = I 
, ol 
Note also that X'x' =I 


In this case x' is said to be orthonormal. 
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APPENDIX 2 


EIGENVALUE PROBLEM WITH A UNSYMMETRIC: GENERAL CASE 
If A is real but unsymmetric its eigenvectors are not orthogonal 
(self-adjoint) as in Appendix 1. However A and *K have the same eigenvalues 
(assumed to be discrete) but different eigenvectors, v; and u,; respectively 


[2]. The n solutions for each may be combined as follows: 


AU = UD Coals 
AV = VD Ce 


where (see Appendix 4 for definition of terms): 


A =n xn real unsymmetric matrix 


U = Modal n x n matrix of eigenvectors u, associated with A 
V = Modal n xn matrix of eigenvectors vj associated with A 
D = Spectral (diagonal) matrix of eigenvalues 


Again as in Appendix 1, after premultiplying (2-1) by V, and also 


after transposing (2-2) and post-multiplying by U: 


{di 
VAU = VUD (2=2) 
and Vau = DW (2-4) 
od 
Putting VU = W and subtracting: 
O = WD - DW 


ie, 


fe. Wij ee 0; 41 # j, since the jy are assumed discrete. Therefore 


Fal 


VU = W (diagonal) (2-5) 
or U and V are said to be biorthogonal. 
Since the basic eigenproblem 
Aui = AGby 
is unchanged by multiplication with a scalar, it follows that any scalar 


multiple of u; is a solution. Hence if U and V are solutions to (2-1) 


LD 





and (2-2), respectively, so also are UE and VF, where E and F are arbitrary 
real diagonal matrices (but note that only post-multiplication by E and F 


is permissible). It is clear that by putting 


BE=Few? (2-6) 


then U'= UE = uw"? 
V'= VF 2 vw? 


will have been mutually normalized, i.e. 
ie > (2-7) 
In this case U' and V' are said to be biorthonormal; they are also said 
to be adjoint. 
It will be noted that if (2-7) is satisfied, then by transposition 
and commutation it is also true that 
uv! = Tv! =v =r 4+ 9'U' (2-8) 
It is also clear that (2-6) is not a unique formula for mutual 


normalization of U and V. The values 
would have done as well, or any of an infinite variety of combinations, 


i.e., there is an infinite multiplicity in the pairs of matrices which 


Satisfy equations (2-8) and which are adjoint eigenvectors of A. 
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APPENDIX 3 


A PARTICULAR EIGENVALUE PROBLEM 


A special case of unSymmetric matrix is derived from the following 
eigenproblem. Let 
AU' = BU'D (3-1) 
where (see Appendix 4 for definition of terms) 


A = real symmetric matrik 


B = real diagonal matrix of non-vanishing diagonal elements 
U' = modal matrix (real) 
D = spectral matrix (real) 


This can be converted to the form of equation (2-1), Appendix 2, since 
yy exists: 
Cu' = U'D (3=2) 
where C = Bes is not symmetric. Therefore, as in Appendix 2, U' is 
not orthonormal, but the adjoint V' exists such that V'U' = I, and 
the property K(C) can be defined as in the text, equations (8). 
On the other hand since A and B in equation (3-1) are both symmetric, 
it is possible by factoring and multiplying through by B 2 to rearrange 


equation (3-1) as follows [4]: 


B aR? B’y = B2U D 
AX xX'* = x'* D (3-3) 
where 
et actin 
A* = symmetric = B “AB 
L 

X'* = orthonormal pseudomode = B7*U' (3-4) 
D = spectral matrix, unchanged 


As in Appendix 1, since A* is symmetric then X'* must be orthonormal: 


ad 
X'*x'* =] 
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It follows that the property N(A*) exists for X'*, as in the text, 
equation (4). But by equation (3-4) this property N(A*) can also be 
associated with U. Since the property K(C) is also associated with U, 
then the question is suggested of how N(A*) and K(C) may be related. This 

% 


question can be amplified as follows. Since B, or B?, or B- » is arbitrary 


(diagonal), then equation (3-4) may as well be written 


uU' = Gx'* (3-5) 
where G =B (3-6) 
= arbitrary diagonal matrix 
The conversion back to the context of equation (3-1), for some given G, 
may be carried out at any time. The general question can now be stated, 
namely, given X'*, an orthonormal matrix with skew-symmetric property 
N(X'*), pre-multiplied as in (3-5) by any G: does the property K belong- 


ing to U' = GX'* have the canonical form sought in the text, for example 


does the main diagonal of K remain empty? 
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Modal matrix . 
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Symmetric........ 


Skew-symmetric... 


Spectral matrix.. 


Orthogonal... 3... 


Orthonormal...... 


Biorthogonal..... 


Biorthonormal.... 


APPENDIX 4 
DEFINITION OF TERMS 


; h 
A quantity in an nf order system represented by a column 
matrix: 


..A composition of all possible eigenvectors into one square 


matrix. The columns of eigenvectors are formed side by 
side. 

ww ; .th 
The transpose of a matrix A, written A, has for its i 
row the ith column of A and for its jth column the jth 
row of A, for all rows i and columns j. ‘4 is said to be 
obtained from A by interchanging rows and columns. 


If A = A, A is said to be symmetric. A symmetric matrix 
must be square and its elements satisfy 4ij = 45> for 
all iand j. The matrix is symmetric with respect to 
its main diagonal. 


, aed 
If the square matrix A is such that A = -A, it is said 
to be skew-symmetric. In a skew-symmetric matrix, 


= -a,,. In particular, the elements a;,; of the 


as: i 


aaa diagonal must be Zero. 


A diagonal matrix with eigenvalues ()\;) as its main 
diagonal. 


If the product of a real matrix X with its transpose is 
a diagonal matrix, W, then X is said to be orthogonal: 
XX = W. Each vector (column) of X is orthogonal to ail 
other vectors (columns) of X. 


A non-singular matrix X is said to be orthonormal if 
a = X; thus an orthogonal matrix X which has been 


normalized is said to be orthonormal: X%X = I. 
U and V are biorthogonal if VU =W (diagonal), 


Lt W = I, then V and U are biorthonormal. 
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APPENDIX 5 
PROGRAM REVEIG 
The first program was based on the analysis of Appendix 2, to solve 


for U and V from a given A: 


AU = UD (5-1) 
AV = VD (5-2) 
where A = Square matrix up to order 25x25. 
U = Modal matrix of eigenvectors with respect to A. 
V = Modal matrix of eigenvectors with respect to A. 
D = Diagonal matrix of characteristic eigenvalues. 


The subroutine used for the eigenvalues is the standard EIG3 subroutine, 
which solves for the roots of an aoe degree polynomial. The eigenvectors 
are not supplied by thisb standard subroutine. It was pecessary to program 
for them a solution of the homogeneous equations A = )\I = 0 for each 
eigenvalue; this was facilited by arbitrarily setting one coordinate of 
each eigenvector equal to unity. The program then tests each diagonal 
element of W = WU for algebraic sign. If W contains negative elements, 
the sign of this element and of the respective column of U is reversed, 
before U is printed out, to permit real square roots of W in accordance 
with equation (2-6) of Appendix 2, in the next step. The change of sign 
is permissible because only direction of the eigenvector is affected 
thereby. 
Biorthonormality was then effected by determining a new U and \V, 

called UP and VP in the print-out, as follows: 

(UP) = UE (5-3) 

(VP) = VF (5-4) 


-k 
where E and F both equal W *. The product VU was then printed out using 


20 





the new UP and VP to check for biorthonormality. 
Following equations (8) of the text the K matrix and its transpose 
were then determined from U' and V', respectively, using the following 


relationships 


K = (rtup) }(1-UP) (5-5) 


ie (vett) (vP-1) (5-6) 


These must, of course, be consistent derivations of K; it was partially 
the purpose of this investigation to verify this fact numerically. 

The program was tested with both a 3x3 and a 4x4 symmetric A, in 
which case UP = VP as was expected, in accordance with equations (9) of 


the text. 


21 





Me ELEN 90 22000 
TMsL oF 
PROGRAM REVEIG 
ODIMENSION A(25325) 3U(25925) sAT(25325) sV( 259 25) sRTR(25) sRT1(25) sAUX 
ee 829) 26 (25 Vek Cea UP (25520 eR 2 50 2 ee) ee 
225) »9FKT(25925) sAUX1 (25925) sRHS(25) 7 
EQUIVALENCE (UP ,W) : \ 
GIVEN MATRIX As FIND THE U MATRIX(ASSOCIATED EIGENVECTORS) AND D(DIAGQNAL MAT 
IX OF ASSOCIATED EIGENVALUES) USING EIG3 SUBROUTINE WHICH BASICALLY USES NTH 
EGREE POLYNOMIALS TO SOLVE THE PROBLEMe EIG3 ALSO RETAINS THE REAL ROOTS AND 
[IMAGINARY ROOTS FOR LATER USE IF REQUIREDe 
} READ 999—5N 
999 FORMAT (15) 
3 Beet N) 23253 
met OP 
3 DO10O IT=ls5N 
READ 998s5(A(I5sJ) sJ=Hl5N) 
10 PRINT9O98Bs(A( 1 9J) sJ=Hl oN) 
998 FORMAT (4E20.210) 
DO 15 I=l>s5N 
DO 15 J=l »N 
15 AUX(IsJ)=A(TsJ) 
CALL EIG3 (AUXsNsNoRTRsRTI 525) 
EIG3 WILL PRINT INTERMEDIATE RESULTS AND RTR>5RTI ° 
DO 11 I=lo5N 
BAD SF(RTI(CT))-LeE-S)lllsllsl2 
MeeerRINT' 997 
997 FORMAT(29HOEIG3 SAYS EIGVAL IS COMPLEX.) 
| STOP. 
Pie 1 (1)=0- 
11 CONTINUE 
REARRANGE EIGENVALUES SO THAs SMALLEST APPEARS FIRST 
Reem MATRIX SINCE EIG3 SOLVES FOR EIGENVALUES ONLY- 
ASSUME ONE COORDINATE OF EACH SET OF EIGENVECTORS EQUALS ONEe: USE EQUATION A 
LAMBDA ¥*I=0 
CALL RFRTR(NsRTR) 
Boei3 I[=15N 
DO 14 J=loN 
DO 14 K=l1,5N 
AUX(JsK)=Al( JU 5K) 
IF (J-K) 145130914 “ 
L30 AUX(JsK)=AUX(JsKI-RTR(IS .- 
~14 CONTINUE 
NM1=N—1 
DO 140 K=1l1>5NMl 
40 RHS(K)=-AUX(KsN) 
! CALL MATALG( AUX »RHS sNM1 9190sDET 925) 
16 DO 18 J=1 5NM1 
emo tJs1)=RHS( J) 
eiNslij=l. 
13 CONTINUE 
PRINT 995 
195 FORMAT (9HOU MATRIX/) 
DO 150 I=l1s5N 
50 PRINT 99lsIs(UlIsJ)sJ=15N) 
CALL MATRA (NsNsAsAT 9259825) 











va 
Re 


Jags 


= ee ee ae —— 





00 19 I=19N 
DO 19 J=loN 
MOP AUX(IsJ)=AT(I J) 
: DO SAME AS BEFORE USING AT VICE As VVICE U c i 
ALSO TRANSPOSE A AND DETERMINE THE V MATRIX OF EIGENVECTORS AND THE ee DIA 
"NAL MATRIX OF EIGENVALUES USING SAME E1G3 SUBROUTINE. 
CALL EIG3(AUXsNoNsRTR9RT1 925) 
DO 20 I=lsN 
| IF (ABSF (RTI (1))-10E-5)120920921 | 
21 PRINT 994 
994 FORMAT(29HOEIG3 SAYS EIGVAL IS COMPLEXe) 
STOP . ) 
Mee RTI(1)=0. 
20 CONTINUE 
REARRANGE EIGENVALUES SO THAT SMALLEST APPEARS FIRST 
CALL RERTR(NsRTR) 
DO 22 J=I15N 
DO 23 J=1>sN 
DO 23 K=l>5N 
BPOXtIsK)=AT(IsK} 
IF (JrK) 239230923. 
23C AUX(JsK)=AUX(JsK)“RTR(1) 
23 CONTINUE 
DO 231 K=1 5NM1 
231 RHS(K)=-AUX(K9N) 
CALL MATALG(AUXsRHS3NM191s0sDET925) 
25 DO 27 J=1lsNM1 
fee ( Js 1) =RHS( J) 
V(NeT)=Hle 
22 CONTINUE 
PRINT 992 
992 FORMAT(SHOV MATRIX/) 
Bee220 I=1lsN 
220 PRINT 99lsIs(Vils Fieoe 1sN) 
_ TEST v TRANSPOSE TIMES U TO SEE IF OFF-DIAGONALS ARE ZERO WHICH ENSURES ORTH 


. ONALITYe CALL RESULT We 


CALL MATRA(NsNsV9AUX 925925) 
96 CALL MAMUL (NsNoNsAUX sU9W925925925) 
PRINT 990 
990 FORMAT (9HOW MATRIX/) 
DO 28 I=15N ; 
28 PRINT 99lsIs(W(IsJ)sJ=15N) 
991 FORMAT(15/(1X7E17e10)) ° 
ASSUME W IS A DIAGONAL MATRIX 
. ARRANGE W SO THAT THE DIAGONALS ARE ALWAYS POSITIVE BY CHANGING SIGNS OF RES 
. CTIVE COLUMNS OF THE U MATRIXe PRINT U ANX W MATRICES AGAINe 
LLL=1 ' = 
DO 50 L=15N 
IF (W(LsL)) 515950950 
Me DO 52 LL=15N 
92 U(LLslL) =-U( LLL) 
-LLL=2 
50 CONTINUE 
GO TO (54,53) >.LLL 
593 PRINT 995 
DO 55 J=l1ls5N 
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Bo PRINT, 99lsIs(UCTsJ)sJ=15N) 
GO TO 56 
C TAKE MINUS SQUARE ROOT OF W AND STORE IN E AND Fe 
154 DO 30 I=l1sN 
. WW=W( Is) 
EF = SIGNF (SQRTF(ABSF (WW) ) sWW) 


EF=1./FF 
E(1)=EF 
30 F(I)=FF 


G POSTMULTIPLY U BY E AND V BY F AND PRINT U*E AS UP AND V¥F aS VP e 
‘ CALL MAMUL1(NsNoUsE 9UP 925925) 
PRINT 980 
980 FORMAT (15HOU-PRIME MATRIX/) 
DO 31 I=l15N 
31 PRINT 991 5I1s(UP(T9J) sJ=l5N) 
CALL MAMULI(NoeNoVeFsVP525925) 
PRINT 98] 
981 FORMAT (15HOV=PRIME MATRIX/) 
DO 32 I=1sN 
32 PRINT 991losIs(VP(I9J) sJ=lsN) 
C TRANSPOSE VP | : 
| CALL MATRA(NoNo VP 9 VPT 925925) 
Cee erIPLY VPT TIMES UP TO CHECK FOR ORTHONORMALITY BY SEEING IF RESULT IS. IDES 
CITY MATRIX. 
CALL MAMUL(NoNo No VPT 9 UP 5 AUX 925525925) 
PRINT 982 
982 FORMAT (322HOV-PRIME-TRANSPOSE TIMES U- PRIME/) 
C COMPUTE K, AND KT MATRIX WITH GIVEN UP AND VP USING THE FOLLOWING EQUATIONS» « 
C (I1+UP) INVERSE TIMFS (I-UP) AND KT= (VP+I) INVERSE TIMES (VP-1) 
wo 33 [=1>5N 
33 PRINT 991leoI9(AUX(I 9J) sJ=l19N) 
DO 40 I=ls5N 
DO 40 J=loN 
IF (I-J) 4139414, 413 
414 AT(IsJ)=let+UP(I 9J) 
BO 10 40 
mis AT(1sJ)=+UP( IJ) 
40 CONTINUE on 
C USE AT AS TEMPORARY STORAGE 
CALL MATALG(AT sAUX1 sNoNo1sDET 25) 
DO 41 I=1>5N 
| DO 41 J=loN 
41 AUX(I,J)= -UP(IsJ) 
DO 42 I=1sN 
42 AUX(I91) =AUX(I51)+1le 
CALL MAMUL(NoNoNsAUXI1 9 AUX sFK 925925925) 
PRINT 983 
983 FORMAT (9HOK MATRIX/) 
DO 43 IT=15N 
43 PRINT 991 sIs(FK(I9J) sJ=l oN) 
, DO 44 TI=l15N 
= DO 44 J=lsN 
IF (I-J) 4155941659415 
416 AT(IsJ)=let+VP(I9J) 
GO TO 44 
415 AT(IsJ)=VP(I9J) 
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CONTINUE 
AS TEMPORARY STORAGE 


CALL MATALG(AT »sAUX1 »NoNo19DET 925) 
DO 45 I=1,N 
DO 45 J=l>s5N 
AUX (IsJ)=VP( TI 9J) 
DO 46 I=15N 
AUX( Is I )=AUX( Ts T)—-le | 
CALL MAMUL(NoNagNoAUX1]-sAUX 9 FKT 925925925) 
PRINT 984 
FORMAT (1LOHOKT MATRIX/) 
DO 47 I=15N 
PRINT 991loIs(FKT( I 9J) sJ=l9N) 
moO 10 1 
END 
SUBROUTINE EIG3(AsNoMsRTRoRTI 9 NQ) 
474 10152 EIGENVALUES OF REAL MATRIGes 


Meee EIGENVALUES OF NON-SYMMETRIC MATRICES 


10 


ia 


eal 


DIMENSIONA(NQ 9NQ) »NC(100) sRTR(M) »RTI(M) 
CALL OVFSET 

TRACE=A(151) 

me 10 I[=2.N 

TRACE=TRACE+A(I 91) 

WRITE OUTPUT TAPE 694 sTRACE 

CALL TRING(AsleE-7T9NoNC» NQ) 

TRACE=A (191) 

BO 11l.1=25N 

TRACE=TRACE+A(1 91) 

WRITE OUTPUT TAPE 694 sTRACE 

NU=0 

NV=0 

IF (NV-N)14912914 
NV=NV+1 

NU=NV | 

IF (NC(NV))15917915 
NV=ENV4+1 

GO TO 16 

IF (NV-NU) 19918919 
RTR(NU) =A(NUsNU) 


f 


RTI(NU)=0. 

WRITE OUTPUT TAPE TENN TS RESIS 
fo TO 13 

IF (NV-=NU-1)20921920 


NP=XMINOF (MyNV) 

CALL LAGER (As1¢E-4sNP»NUsNVoRTR RTI 9 NQ) 
Beto 13 

RR= 05% (A(NUsNU)+A( NV NV) ) 

EL=RR¥*2—A (NU sNU) ¥A(NV NV) +A(NUSNV) ¥ACNV 9 NU) 
S=SQRTF (ABSF(E1)) 

MEE 1)22923523 

RTR(NU)=RR+S 

RTI (NU)=0. 

RTR(NV) =RR-S 

RTI (NV) =0. 


x 


WRITE OUTPUT TAPE 695 sRTR(NU) RTI (NU) #RTR(NV) 9RTI (NV) 


GO PO. 13 25 


j 
‘ 


FIG300C 
E1G300( 


ADDED 
E1G300( 
E1G300( 
E1G300( 


“£16500 


E1G300: 
E1G300: 
E1G300: 
E1G300: 
E1G300: 
E1G300: 
E1G300: 
E1G300: 
E1G300: 
E1G300: 
E1G300z 
E1G300z 
E1G300: 
E1G300: 
E1G300; 
E1G300: 
E1G300: 
EI1G300; 
E1G300: 


E1G300% 
E1G300: 

E1G300: 

E1G300: 
E1G300: 
E1G300: 
E1G300: 
E1G3002 
E1G300% 
E1G300: 
E1G300¢ 





Be RTR(NU)=RR 


RTI (NU) =S 
RTR(NV)=RR 
RTI (NV) =<S 
GO TO 25 
2 X=Oe 
G CALL FPOLD 


DO 24 J=l»M : 
2¢ X=X+RTR(J) 
WRITE OUTPUT TAPE 656 9X 
© §6RETURN 
 & FORMAT(1HO48X»7HTRACE =E1608) 
' 5 FORMAT (11HOEIGENVALUE 12X»2E2008) 
é FORMAT(1HO35X920HSUM OF EIGENVALUES =E1608) 
END 
CO51/2ALMOST TRIANGULAR (HESSENBERG) SUBROUTINE 
SUBROUTINE TRING(AsEPSsNoINT» NQ) 
DIMENSION A(NQsNQ) 9 INTCNQ) 
WRITE OUTPUT TAPE 691 
N1=N-1 
N2=N-2 
00 21 J=1sN1 
S=ABSF(A(JsJ+1)) 
mi =J+] 
J2=J+2 
1 
NJ1=N-J1 
IF (NJ1) 1551556 oe, 
6 DO 12 K=J25N 
T=ABSF(A(JUsK)) 
IF (T-S)12912911 
ial eK 
S=T 
12. (CONTINUE 
. IF (L-J1)13515913 
13 MO 131 K=lsN 
TH=A(KyJ41) 
Bt Ky J+ 1)=A(KoL) 
Be At KoL)=T 
14: BO 14] K=15N 
‘THA (J+1 9K) 
AC Jt+1s K)=A(L oK) 
141 A(LsK)=T 
15 IF (S-EPS*MINIF(ABSF(A(J9J) ) sABSF(A(J+159J+1))))16916917 
16 0 
NJ1=0 
GO TO 181 
7 THA JsJ4+1) 
BPO 18 K=J25N 
18 Al JSsK)=Al Js K)/T 
181 'DO 20 I=l »N 
~~ ;M=MINOF (Js 1-2) 
U=O0e6 i 
Mee (NJ1) 1951937 ; 
7 DO 8 K=J25N 
8 “USUFA(KsI)*¥A(J9K) 


{ 
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‘st 


EG tire 
E1G2Ce 
B1GSOR 
EI1G70G 


~ EIGS0C 
EG 3 


EIG30C 
E1G2o0 
E1G300 
ELG20G 
EIG300 
EIG300 
ELGSc6 
E1GI0¢ 
ElG2or 
ELG20e 


E1Gs0C 
EIG30¢C 
ElGaQ6 
ETG30C 
EIGICG 
EI1GS06 
ElG2Oe 
EIG206 
E1G20G 
ElG30¢ 
EIG30¢ 
EITG30C 
EIG30C 
ETG30¢ 
E1G206 
ElG@Ce 
EIG30C 
ELGZO¢ 
EIG300 
EI1G30€ 
ETG3CEé 
E1Ge20¢ 
EIG306 
EIG30C 
EIG306 
EI1G36C 
ELG3 OC 
E1GZCe 
EIG3GC 
EIG30C 
EIG300 
E1G20C 
ENGSCC 
E1G3C0 
EIGICC 
E1G300 
EGS CC 
ELGS GG 





IF (M) 2052099 

DO 10 K=1>5M 

USU-A( Kol) *A(J4+19K41) 

A( J+ls1)=A(J+1l91)+4U 

INT( J) EL 

BT iN)=0 

RETURN 

FORMAT ( 1HO48X » 22HALMOST TRIANGULAR FORM)’ 
END - 


061/3LAGUERRE METHOD 


‘ 
‘ 


\ 


13 


4 


5 


i? 


SUBROUTINE LAGER(AsEPS  sN1l»NUsNoRTRoRTI »NDI ) 
DIMENSION A(NDI »NDI) 9P (69101) sRTR(NDI) sRTI(NDI) 98(6) 
WRITE OUTPUT TAPE 691 

ONCE=0.6 

mela le 

NUQ=NU-1 | 

iy =0 

DELOLD=1le 

ROLD=1. 

EGSUM1=06 

EGSUM2=0-6 . 
CALL FPTEST(Z) 

DO 10 L=256 

ote NU)=0. 

NU1=NU+1 

CUP=0. 

DO 11 J=NU1>5N 

CUP=CUP+ABSF(A(J—-1l9J)) 

CUP =CUP/FLOATF(N-NU) 

mAP=0. 

P(1sNU)=1. 


FIND TRACE OF H AND H SQUARED 


SPUR1=A (NU» NU) 
SPUR2=A(NU»sNU) ¥¥2 

DO 13 J=NU1oN 

SPURI=SPUR1+A(JoJ) 

SPUR2=SPUR2+A(J9J) ¥¥24+20¥*A(J—19J) ¥A(J9J=1) 


INITIAL ITERATE FROM INFINITY 


eek=tGSUM1-SPUR]1 
S2R=SPUR2-EGSUM2 
F1l=N-NUQ 

IF (ABSF (S1R)+ABSF(S2R)-1- E-7#CAP)15, 155 16° 
XBAR=CUP 

YBAR=0-6 

io 7O 23 

F2=Fl-le 
DR=F2*(F1*¥S2R-S1R¥*2 ) 
ER=SQRTF(ABSF (DR) ) 

ITF (DR)17918518 
XBAR=-2e*S1R/F 1 
‘YBAR=2.*ER/F1 


m0 TO 23 
| Zi 


- E1G368 


E1G3¢@ 
E1G3¢@ 


ate 16307 


FEIG300 


7 


EIG308 
E1GaGal 


£1G301 


ELG2Z0u: 
EIG2Ga 


EIG301. 
EIG208 
EIG301 
EIG301- 
ETG301 
EIG301. 
EIG301 
EIG301 
ETPGe Ge 
oe 
EIG30@7 
EPGsen" 
ElGooue 
ETGS0n3 
E1[G3609 
ETGlo1 
EIG301. 
E1G301. 
EIG30l; 
EI1G305 
EIG30 18 
EIG3049 
EIG308 


E1G301? 


E1G301. 
E1G30 18 
EI1G301: 
E1G301: 
EIG301: 
E1G301: 
EIG301: 
EI1G301< 
E1G301: 
E1G301¢ 
E1G301¢ 
E1G301: 
EIG301¢ 
E1G301¢ 
EIG301¢ 
E1G3014 
E1G301¢ 
EIG3015 
EIG301! 
EIG301! 





18 Y3SAR=0. 
F2=SIGNF(1.5S1R) 
IF (S1R)22s20922 
mZ20 F2=064 - P 
= 22 XBAR=-(SIR+F2*ER) /F1 
G 


C 


C EVALUATE POLYNOMIAL AND DERIVATIVES 
23 IF (ABSF (YBAR)-ABSF(XBAR)¥*1 66-6) 24525925 
24 YBAR=0. 
25 M=6 
IF (YBAR)27926927 
26 M=3 


27 DO 34 K=NUsN 
T=-A(KsK+1) 

0O 34 L=l»M 
S=SIGNF (le9305-FLOATF(L)) 
M1=L+3*XFIXF(S) : 
R=-XBAR¥P(L9K)+YBAR®S*P (M1 9K) -FLOATF(XMODF (L—-193))*P(L—1 5K) 
BO 28 J=NUsK 

28 R=R+P(LoJ)*#A( Ko J) 
CALL OVFTST (Z) 
1F(Z)29932929 

moe «6Z=06 
P(1sNU)=1.E-10*P( 1 5NU) 
TF (P(1lsNU)) 30230927 

30 F=FLOATF(K—NU) /FLOATF (N-NU+]1 ) 
WRITE OUTPUT TAPE 6929XBAR»YBAR »P(1 NU) 9F »ONCE 
XBAR=XBAR*¥F 
YBAR=YBAR*F 
Bw ,NU)=1.0 


GO TO 27 
32. IF (N-K) 33533934 : 
a 1=1.0 


34 P(L»sK+1)=R/T 
SCALE DOWN : 


DO 39 K=1s6 
39 B(K)=06 
DO 35 J=l»M 
peo t(J)=P(J»sN+]1) 
Gl= ABSF(B(1))+ABSF(B(4)). 
'G2= ABSF(B(2))+ABSF(B(5)) 
G3=ABSF (B(3))+ABSF(B(6)) 
WRITE OUTPUT TAPE 6929XBAR sYBAR sG1 9G29G3 
Z FORMAT (8H ITERATE20X9E1508 95XsE150898X93E1504) 
D=ABSF(B(1)) 
DO 36 K=29M 
36 D=MAXI1F(DsABSF(B(K))) 
(CALL SCALE(D»sB 9M) 
air (G1) 41541543 
l 


‘REMOVE KNOWN ROOTS 


C 28 
43 Q1R=06 


et 8, SS NR Et de. ptthinliny TN ee tet ven AOS tN eed Sa Ait amet mY + tog me —* are 


EGG) 
E1G3GH 
EIG3Q) 


_ ,£1G301 


EIG301 
ElG3Ga 
EIG3@] 
EIG3G1 
EIG3@) 
EI1G3@i 
EIG30] 
E[G2@I 
EIG2@1 
EI1G3G) 
EIG30) 
EIG3G1 
E]G2@l 
EIG3@) 
EtGacs 
EIGI¢i 
E1G3@) 
ADDED 


-ELGS@% 


ELGeG 
EIG303 
EIG36@ i 
EIGZG} 
EI1G301 
EI1G3@1 
E1GaG) 
ElGle! 
E1lGs@} 
EI1G30) 
EIG3¢% 
EIG30] 
EIG2¢e) 
E1G3q1 
E1G3e3 
ElG365 
EIG308 
EIG30% 
EIG304 
EIG3G% 
EIG301 
EIG30] 
EI1G309 
EIG303 
EIG301 
EIG2@4 
EIG304 
EI1G304 
E1G30¢ 
E1GI04a 
EIG304 
EIG302 


—IG304 








Zl 


44 


19 


45 


46 


i) 


a2 


Gal I =O, 

Q2R=0. 

e21=0- ; 
IF (NUQ=NU) 19921921 
DO 44 J=NUsNUQ 
DI=RTR(U)=XBAR 
D2=RTI(J)-YBAR 
D=D]1*#24+D2* #2 
D1l=D1/0 

D2=-D2/D 
Q1LR=Q1R+D1 

Ol =Q11+D2 
Q2R=Q2R+D01**2-D2¥**2 
Q21=Q21+2-.*D1*D2 


FIND S1 AND S2 


TIR=B(2)/B3(1) 

T11=0. 

f2R=B(3)/38( 2) 

mZ21=0. 

IF (YBAR)45346945 

DI=B(1)**2+B(4) #¥2 
D2=B(2)**2+B(5) *#*2 
T1IR=(B(2)*B(1)+B(5)*B(4))/01 

je =(8(5)*6(1)-B(4)*B(2))/01 
T2R=(B(3)*B(2)4B(6)*B(5))/D2 
T21=(B8(6)*B(2)-B(5)*B(3))/D2 
S$1R=T1IR+Q1R 

SlI=HT1LI+Q11 
Pere=TIR*¥(TIR=-T2R)-TLIVM(T1II-T21)-Q2R 
S2T=TIR¥(TLI-T21)4+T11*(TIR-T2R)-Q2 1 


menD THE NEXT ITERATE 


my =LLY+1 
D=ABSF (XBAR)+ABSF(YBAR) 


IF (1¢E+7-D* (ABSF(S1R)+ABSF(S11I))) 41941942 | 


MARK=] 

GO TO 100 

G=N-NUQ . 
IF (YBAR-ABSF(X))503950949 - 
S1LI=SlI14+16/(2e¢*YBAR) 
S2R=S2R+10/(4e*YBAR¥X2 ) 
G=G-l. | 

IF (BL1)65965 966 

H= 0 5#(G=2.6) 

GO TO 67 

H=G-le 
DR=H*(G*S2R-S1R¥¥24+S511%**2) 
DI=H*(G¥S21-26*S1R*S11) 

IF (DI)53951553 

EI=0-6 

ER=SQRTF(ABSF(DR) ) 

IF (0R)52954,554 

EI=ER 


29 


| 
E1G302C 
E1G3020 
E1G3021 


EIG@GZe 


E1G3021 
F1G3021 


“ELGAQ2Ar 


E1G302] 
E1G302) 
ETG362) 
ETG3072 i 
F1G3021 
ElG3022 


E1G3022 


E1G3022 
EIG302z 
EIG3022 
E1G3022 
E1G3022 
E1G3022 
E1G3022 
E1G3022 
E1G302? 
E1G30223 
E1G30223 
E1G3023 
E1G3023 
E1G3023 
E1G3023 
E1G3023 
E1G3023 
E1G3023 
E1G3024 
E1G3024 
E1G3024 
E1G3024 
E1G3024 
EIG3024 
E1G3024 
E1G3024 
EIG3024 


EIG3024 


E1G3025 
EIG362Z5 
E1G3075 
E1G3022 
EIl1G3025 
EPG 302 
ElLG3025 
BUG Sie) 25 
EI1G3025 
ELGeoze 
E1G63026 
E1G30:26 
EIG302¢ 
ELGSoze 
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7a 
f 2 
n3 





FR=0,. 

GO TO 54 

MALL CXSQRT(BResDISERsSET) 
ME (SIRKER+SL IE!) 555956956 
ER=-ER 

eI =-E! 

OL=S1IR+ER 

D2=SLI+E! 

D=D1**24+D2**2 

X=-G*D1/0D 

XBAR=XBAR+X 

Y=G*D2/D Y 
YBAR=YBAR+Y 
DELNEW=A8SF(X)+ABSF(Y) 
RNEW=DELNEW/DELOLD 

D=ABSF (XBAR)+ABSF(YBAR) 


TEST FOR LINEAR CONVERGENCE 


Mee el Y-3)62362957 . 

mee (DELNEW=-MAX1F(36*DELOLDs e5*D) 157155719570 
MeetS5L1) 571595713572 

DELOLD=CAP | 

ROLD=3-6 2 
IF (LLY~-15) 1459145100 

IF(RNEW-e 7*ROLD) 62958558 

MARK=3 

IF (DELNEW-eO01*EPS*¥CAP) 70959559 
IF (BL1)61s61560 
XBAR=XBAR-X 

YBAR=YBAR=-Y 

‘BL1=0. 

GO TO 48 

BLi=l. 

m0 TO 63 


i < 
TEST FOR AN EIGENVALUE 


IF (DELNEW-EPS¥MAX1F(D» e001¥CAP) 164564963 
DELOLD=DELNEW 

ROLD=RNEW 

Mette LY-15) 233,233:100 


DO WE HAVE A COMPLEX APPROACH TO A REAL ROOT 


MARK=2 

BL1=1 

IF (YBAR) 713100571 

IF (G2*ABSF (YBAR)—-G1) 7251005100 
mr CONCE ) 735739100 

X=06 

ONCE=1. 

1 YBAR=0. 

aso 10 63 


WE ACCEPT (XBARsYBAR) AS A ROOT 
: 30 


E1G302 
F1G302 
E1G302 
E1G302 
E1G302 
E1G302 
E1G302 
E1G302 
FiGeo2 
E1G302 
E1G302 
E1G30? 
E1G302 
E1G302 
E1G302 
E1G302 
E1G302 
E1G302 
E1G302 
E1G302 
E1G302 
E1G302 
E1G302 
E1G302 
E1G302 
E1G302 
E1G302 
E1G302 
E1G302 


"ETG202 


E1G302 
E1G302 
E1G302 
E1G302 
E1G302 
E1G302 
E1G303 
E1G303 
E1G303 
E1G303 
E1G303 
E1G303 
E1G303 
E1G303 
E1G303 
E1G302 
E1G303 
£16303 
F1G303 
E1G303 
E1G302 
E1G303 
E1G303 
E1G303 
E1G303 
E1G203 








80 
8 


83 
30) 


—6B4 


86 


101 
1 





NUQ=NUQ+] 

RTR(NUQ)=XBAR 

IF (ABSF(YBAR)=~2.001*ABSF(XBAR) ) 74574575 

YBAR=0.6 

IF(NUQ=NU) 99769 9 

mee RTI (NUQ-1)) 76576577 

YBAR=ABSF(YBAR) 

RTI (NUQ)=YBAR > 

CGO 10.78 

RTI (NUQ)=-ABSF(YBAR) 

WRITE OUTPUT TAPE 693sRTR(NUQ) »RTI(NUQ) » LLY» MARK | 

FORMAT (1 LHOEIGENVALUE1 2X 92E2068912X%313917H ITERATIONSsTEST I1//) 
LLY =0 . ° 
CAP=MAX1F(D»sCAP) 

ZEeLOLD=1. 

ROLD=1. 

EGSUM1=EGSUM1+RTR(NUQ) 
EGSUM2=EGSUM2+RTR(NUQ) *#2-RTI(NUQ) ¥¥2 
IF (NUQ-N1)80+1019101 

IF(YBAR)83384,81 

YBAR=-YBAR 

Bo TO 23 : 

IF (NUQ=NU ) 31584531 


—-s 


RTI (NUQ=1)=.5#(RTI(NUQ=1)-RTI(NUQ) } 

RTI (NUQ)==RTI(NUQ-1) 

A NEWTON ITERATE TOWARDS NEXT ROOT a 
ONCE=0. 

Z=0-6 


IF ( (ABSF(Q1R)+ABSF(Q11))*D-10000. ) 85585514 

IF (ABSF (EGSUM1—SPUR1 ) +ABSF (EGSUM2—SPUR2 )-1eE-5*CAP) 15915586 
DR=B(3)+2.*(B(2)*Q1R-B(5)¥*Q11) 
DI=B(6)+2.*(8(2)*Q11+B(5)*Q1R) 

D2=DR¥#¥#24+D1 **2 

XBAR=XBAR-2.*(DR*¥B(2)+DI*B(5))/D2 

YBAR=ABSF (YBAR—26#(DR*B(5)=-DI*B(2))/D2} 

GO TO 23 

RETURN 


EEG Aes 
ElG3oe7 
B1G3034¢ 


“me 1G3032 


E620 3 
EhGa6 32 
E NG eiG 3 


BeBe 219 6) 


E1G303z 
F1G3032 
E1G3032 
E1G3032 
F1G3032 
E1G3032 
E1G3037 
E1G303: 
E1G303: 
E1G3032 
E1G303: 
E1G303: 
E1G3034 
E1G3034 
E1G3034 
E1G3034 
FE 163034 
E 163034 
E1G3034 
E1G3034 
E1G3034 
E1G3034 
E1G3035 
F1G3035 
EIG3035 
E1G3035 
E1G3035 
£1G3035 
E1G3035 
E1G3035 
E1G3035 
E1G203« 


FORMAT ( 1HO050Xs 19HLAGUERRE ITERATIONS/31Xs9HREAL PART1OXs1OHIMAGe PFIG303¢ 


LART22X51HP11Xs7HP PRIME6Xe11HP DBL PRIME) 
END P 


)1/5COMPLEX SQUARE ROOT 


SUBROUTINE CXSQRT(AsBsXsY) 
F=MAX1F (ABSF(A) sABSF(B)) 
F=F¥SQRTF( (A/F) ¥¥2+(B/F)¥*2) * 
Mm (Al1,152 : 
Y=SQRTF ((F-A)*.5) 

X=.5*B/Y 

1F(X)45353 

X=-X 

--Y 

BO TO 3 

X=SQRTF((F+A)%.5) au 
Y=.5*B/X 


ET Ge Oce 
ELGZOaG 
ElG3036 
ElLG3026 
EIG202e 
BING BNO Sie 
EIG3036 
El1G302e 
ElGe Ose 
E1G2037 
E rG3@ray 
E1G3087 
ElGs0e7 
ElGgGe 7 
E1G3027 





3, RETURN FilGso2 
END E1G302 
SUBROUTINE SCALE (0D 5B 5M) : 

DIMENSION B(6) 
S = 2. ** (XINTF (LOGF(D)/-69314718056) +1) / 


DO 1 IT=loM 
Mmea(l) = BII)/S 
RETURN 
END 
IDENT OVFSET 
ENTRY OVFSETsOVFTST »sOVFLU 
@erseT SLJ x 
NOP 
CALL SELECT* 
| al OVFLU 
me CALL ERROR* 
ENA 0 . 
STA =SZIG = 
SLU OVFSET 
OVFLU SLU %% 
EDA =Dle re 
STA Zl 
CALL SELECT* 
1l OVFLU 
CALL ERROR* 
| Sap OVFLU 
OVFTST SLJ %% 
_ LDA x _ 
ARS 24 
SAL % +2 
INA l 
SAU OVFTST | 
m LDA Z1G 
STA %% 
ENA 0 
STA Z1G 
7 SLU OVFTST 
END i , 
SUBROUTINE MATALG(AsXsNR9NV9IDOsDET »NACT) 0000 
DIMENSION A(NACT sNACT) sX(NACT sNACT) 0900 
mre{1DO) ls2ol A : 0000 
1 DO 3 IT=l»sNR a 0000 
wo 4 J=leNR - aie.0 
4 X(15J)=020 0000 
xi 1>1)=1.0 0000 
NV=ENR | = atorale 
EeDET=1.0 ce 
NR1=NR-1 0000 
moe5 K=1,NR1 000] 
IRL=K+1l . ; 0001 
;PIVOT=0-0 0001 
-Z=ABSF(A(IsK)) 0001 
IF(Z-PIVOT) 69697 0001 
PEPIVOT=Z 9901 
IPR=I : 0001 


Zz 





— 


CONTINUE 90018 


IF(PIVOT) 85998 00019 
DET=0.0 0002¢ 
RETURN | 00021 
IF(IPR=K) 10911510 . O0GZ2 
DO 12 J=K»sNR 00022 
Z=A(IPR»J) - 90024 
PIPR,J)=A(KoJ) :; 00025 
mek 9 J) =Z | 00026 
Ba 13 J=1sNV 00027 
m7 =X(IPR»oJ) | 00028 
X(IPR»J)=X(KoJ) 00029 
mk sy J)=Z 00030 
DET=-DET . 00031 
DET=DET¥A(Ko5K) 00032 
PIVOT=1eO/A(KoK) 00033 
DO 14 J=IR1»NR 00034 
Peeks J)=A(Ko J) *¥PIVOT 00035 
DO 14 1=IR1,NR | 00036 
Rls J)=ATsJ)—A( Lo KY *A (Kod) , 00037 
DO 5 J=l»NV : | 00038 
MeiX{(KsJ)) 1555515 00039 
MiKsJ)=X(K»J)*PIVOT 00040 
DO 16 I=IR1l>5NR 90041 
MOTs J) =X( Lo JPmACT 9K) #X (Kod) 00042 
CONTINUE 00043 
IF(A(NRoNR)) 1799917 | 00044 
DET=DET*A(NR»NR) 7 00045 
PIVOT=1e0/A(NRoNR) 00046 
wo 18 J=1>,NV 00047 
X(NR»J)=X(NR9oJ)*PIVOT 00048 
DO 18 K=1 5NR1 00049 
I=NR-K 00050 
SUM=0.0 00051 
DO 19 L=I»NR1 00052 
SUM=SUM+A (I 9 L4+1)*¥X(L+19J) | 00053 
¥(1sJ)=X(1sJ)-SUM ; 00054 


END 00052 
SUBROUTINE MATRA(M Ns Ag B»MD ND) | 


DIMENSION A(MDsND) »B(MD ND) 

DO 3 I=l»M 

Bons J=1.5N 

Sty J)=A(Js1) 

RETURN ~ : 
END 

SUBROUTINE MAMUL(MsNsLb9AsBsC»sMDsNDo9LD) 
DIMENSION A(MD ND) sB(ND LD) s»C(MDyLD) 
DO 2 I=lsM 

wo 2 J= LoL 

S=0-6 

DO 1 K=loN 

S=StA( I 9K) *B( Ks J) 

C(I ,J)=5 

‘RETURN ‘ 

END 


“SUBROUTINE MAMUL1(MsNsAsE2BsMDsND) se 








DIMENSION A(MD»sND) sE(ND) 9B (MD—ND) 
pom i= ly 
DO | J=l1.N 
fs ( | sJ)SACTs J eet) 
RETURN. 
¥ END 

SUBROUTINE RERTRI(NoRTR) 
DIMENSION RTR(25) 
NMl=N=-1 
BO) 12 T=1.NM) 
TP=I+] 
DO 12 J=IPoN 
IF(RTR(IT)-RTR(J)) 12912913 

13 TFMP=RTR(1) 
RTR(I)=RTR(J) 
RTR(J)=TEMP 

12 CONTINUE 
RETURN 
END 

END 
FINIS 
=EXECUTEs3- 
4 
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APPENDIX 6 
PROGRAM REVEIGL 

The second program, REVEIG1, was formulated for a maximum 25x25 
input based upon equations (7) of the text: 

UP = (1-K)(I+k)”! (6-1) 
vP = (14K) (1-%)~* (6-2) 
The input is any real square matrix K. 

In addition, various operations with UP and VP are printed out for 
investigation of any obvious relationships, particularly among the diagonal 
elements; these are Jvi, vv, Tul, u'd', and ¥'u'. 

The K outputs of REVEIG were tested as inputs to REVEIG] in 3x3 and 
4x4 matrices, and in both cases the UP and VP compared with the UP and VP 
in REVEIG. T'u' always resulted in the identity matrix. When K was skew- 


symmetric, UP and VP were equal as was expected. 


35 





pe Reese 22000. 
9 @ 
PROGRAM REVEIGI 
BDIMENSION UP(25525) .VP(255925) sFK(25>5 Zo) eK t2 28 Be TeVET (25925) 
muPT (255925) »sAUK( 25925) »AUX1 (25325) 5T(25525) 
C€ GIVEN A K MATRIX FIND ITS ASSOCIATED U PRIME MATRIX USING THE COE eae EQUA 
@ QN. UP=(I-K)(I+K) INVERSE 
Wee lt READ 8995N 
" IF (N) 29253 
251 OP 
3 PRINT 8705N 
Meo ?O FORMAT (15H1K MATRIX, N=13) 
| 899 FORMAT (15) 
| DO 48 IT=l1>5N 
READ 898s5(FK(I1I59J) sJ=l1 oN) 
48 PRINT 8985(FK(I1I9J) sJ=Hl oN) 
898 FORMAT (4E20.10) 
G STORE (I+K) IN AUX 
DO 49 J[T=loN 
DO 49 J=l>5N 
IF(I-J) 490,491,490 
491 AUX(IsI)=FK(IsI)+1. 
GO TO 49 
490 AUX( 1 sJ)=FK(19J) 
49 CONTINUE 
C USE AUX AS TEMPORARY STORAGE 
MeeNVERSE OF (CI+K) IN AUX1 
CALL MATALG( AUX sAUXI1 »NoNol sDET 225) 
G@eoTORE (I-K) IN T 
DO 50 I=l1>5N 
DO 500 J=15N 
BON i1>+J)=-FK(1sJ) : . 
DO 51 I=loN 
SFC IlsIl)=TtCIsI)+le 
C PREMULTIPLIES (I+K) INVERSE BY (I-K) 
CALL MAMULI(NeNoNoT sAUX1] 9UP 525925525) 
PRINT 897 
897 FORMAT (15HOU=PRIME MATRIX/) ; 
DO 52 IT=1l9N | 
52 PRINT 89691 5(UP(I5J) sJ=loN) 
M26 FORMAT(I5/(1X7E17-19)) 
C TRANSPOSE THE K MATRIX AND FIND ITS ASSOCIATED V MATRIX USING THE FOLLOWING €£ 
M@MUATION. VP=(I+KT)(I-KT) INVERSE | 
| CALL MATRA(NsNogFKoFKT 925925) 
MeotORE .(I+KT) IN AUX 
DO 53 I=l»N 
DO 53 J=l>5N 
AUX(I,J) =FKT(1I5J) 
MeeTORE {I-KT) IN AUX] | 
AUX1(IsJ)=-FKT( I,J) 
Met I-~-J) 535543;53 
~. 54 AUX( 151) = AUX(IsI)+le 
AUX1L( I sT)=le-FKT( I 591) 
593 CONTINUE 
C STORE INVERSE OF (I-KT) IN T 


CALL MATALG(AUX1 »TsNoNol sDET 925) 36 





PREMULTIPLY (I-KT) 


CALL MAMUL(UN sNoNoAUX 5 T 


me INT 895 


INVERSE BY 
9VP 925925925) 


Clit) 


995 FORMAT (1S5HOV-PRIME MATRIX/) 


BS 


PRINT OUT ALSO V PRIME TRANSPOSE TIMES V PRIME, 
SAME FOR U PRIME AND U PRIME TRANSPOSE. 


Eo 

B94 
893 
892 
B91 


CHECK MOO OEE IF V PRIME TRANSPOSE=1 WHICH ENONIES ORTHONORMALITY« 


890 


56 


om | 


58 


Bo 


moe (Z-PlVOT) 


BO 55° 1=1;N 
PRINT 8985 I>» 


(VP( IJ) > 


FORMAT( 9HOVPT X VP / 
FORMAT( SHOVP X VPT / 
FORMAT ( S9HOUPT X UP / 
FORMAT ( 9HOUP X UPT / 


FORMAT (9HOVPT X UP/) 
CALL MATRA (NsNoVP9VPT 925925) 


CALL MATRA (NoNo UP, UPT » 
CALL MAMUL (NosNoNoVPTsVPo5T 325925925) 


PRINT 894 


MmeolGN 56 TO ISW 


GO TO 60 


CALL MAMUL (NsNsNoVPo9VPT 97525925925) 


BeINT 893 


ASSIGN 57 TO ISW 


GO TO 60 


CALL MAMUL (NoNsNoUPT sUP9T 925925525) 


mei NT 892 


ASSIGN 58 TO ISW 


so TO 60 


J=15N) 


) 
) 


) 
) 


25925) 


CALL MAMUL(NoNsNoUPsUPT 57925925525) 


meoht 89) 


ASSIGN 59 TO ISW 


mo 10 60 


CALL MAMUL(N»Ny Ny VPT sUP» T25wies +251) 


PRINT 890 


ASSIGN 71 TO ISW 


DO 61 I=l5N 
PRINT 8965 I+ 
GO TO ISW 

oo TO 1 

END 


SUBROUTINE MATALG(AsXsNRsNV3ID0O5DETsNACT) 


(T(Il sd) sJ=l oN) 


f 


DIMENSION A(NACT 4 NACT) sX(NACT »NACT) 


Ment DO) ls2s1 
DO 3 IT=1>5NR 
DO 4 J=l NR 
X(IsJ)=0-0 
m1 91)=1-0 
NV=NR 

ee(=)1.0 
NR1=NR-1 

DO 5 K=1sNR1 
IRL=K+1 
PITVOT=0.0 

DO 6 IT=KsNR 
ZL=ABSF(A(IsK)) 


69697 


ai 


V PRIME TIMES’ V PRIME: TRANSP 


OOCOE 
0000] 
N000¢ 
0000: 
00004 
00005 
O0O00E 
0000% 
00008 
00005 
OUO1S 
QOO013 
O001< 
OOOGAG 
00014 
GOO 








‘12 


ES 


a1 


r9 
18 


PIVOT=Z 

IPR=] 

CONTINUE 

IF(PIVOT) 859938 

DET=0.0 

RETURN 

me IPR-K) 10514510 

DO 12 J=KsNR . 
m=A(IPRsJ) 

A(IPRsJ)=A(Ko9J) 

A(K,J)=Z 

DO 13 J=elsNV 

Z=X(IPRoJ) 

XC IPRsJ)=X(KoJ) 

X(Ket)=Z 

mET=—-DET 

DET=DET#*¥A(KsK) 
PIVOT=1e0/A(K 9K) 

DO 14 J=IR1>5NR 
A(KoJ)=A(KsJ)*PIVOT 

DO 14 JT=IR1>5NR 
A(ITsJ)=H=ACITsJI-AC IT sK)#A(K 9J) 
DO 5 J=lsNV 

met X{(KoJ)) 15235515 

Merk 9s J)=X(Ke9J)*PIVOT 

DO 16 I=IR1»sNR 
X(LysJV=XC IL sJ)mAC IT 9K) *X (Kod) 
CONTINUE 

IF(CA(NRoNR)) 1799917 
DET=DET#A(NRoNR) 
PITVOT=1eO/A(NRs NR) 

DO 18 J=FlsNV 
X(NR»sJ)=X(NRoJ) PIVOT 

DO 18 K=l»NR1 

I=NR-K 

SUM=0.20 

DO 19 L=I»sNR1 
SUM=SUM+A( I yL+1)*X(L+1 9J) 
X(19J)=X(19J)-SUM 

END 

SUBROUTINE MATRA(M 9N 9A 9B MD 5ND) 
DIMENSION A(MD ND) 9B (MD sND) 
DO 3 :T=1l9M 

DO 3 J=El>sN 

B(IsJ)=Al(JoI) 

RETURN 

END 

SUBROUTINE MAMUL(M»Nol sAsBesCoMDsNDyLD) 
DIMENSION A(MD sND) sB(NDoLD) »C(MDsLD) 
DO 2 I=lsM | 
DO 2 J=losL 
S=0-6 

DO 1 K=loN 
S=StA( I sK)*B(KsJ) 
C(IsJ)=S 
RETURN 


} 
i 
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fd 


*! 


OOO t 
0001: 
SUCH 
renee: 
0002¢ 
0002) 
00022 
O00 2 
00024 
COZ 
DOUZE 
00027 
OO002E 
OOH 8) 2. 
O003G 
00034 
000 3< 
O00 5: 
0003¢ 
DAG IONS, 2 
000 3€ 
00037 
HOO SE 
O00 35 
N004C 
0004] 
0004. 
0004. 
00044 
QOO04E 
OO04E 
00047 
NOO4E 
0004S 
O005C 
00037 
0005¢ 
O18 (6) & 
00054 
Q002e 





END 
END 
FINIS 


See CUTE 9 Bie 
~ & 


¢ 


ee 


A.ANNNNANNNNDOE 
1.QN0900090000E 
4,.NN0QNNN0000F 
SePOOQQQ0000E 


00 
00 
On 
00 


2eDD0DDND0DNN00OEF 


0.00000N000NE. 


ZeINDNNDONHNONE 
4 DODDDNNNNNE 


ole) 
An 
ole 
00 


Si) 


5 ONNNNNONNNE 
2 e00000NN0NNE 
0.0000000N00F 


. ZeONDDNNDNDONNE 


LeNNNNNADNNNDNE 
2 eNNNNNBNNNNE 
3.9nnNONDAANE 
O.aNND0NNNNNE 





APPENDIX /7 
PROGRAM REVEIG2 
A short third program, REVEIG2, was formulated to carry out the cal- 

culation in equation (3-5) of Appendix 3, and to compute the corresponding 
matrix K(U') from equation (8) of the text. The program printout notation 
requires special description however. The input to REVEIG2 was an arbitrary 
real diagonal matrix, G, and an orthonormal matrix X'*, as in equation (3-5). 
The latter matrix was, however, the result obtained from REVEIG for a sym- 
metric matrix, and called U' or UP in REVEIG instead of X'*. It is there- 
fore called UP in REVEIG2, also, necessitating the introduction of the 
program description UPP, or U'', for the calculation in equation (3-5). 


Thus in REVEIG2 the following is calculated 


G(UP) (7-1) 


URE 


and 


~ 
li 


(1+UPP) + (1-UPP) (7-2) 


40 





meOOP s,s 5YOST 9S/15/2533,50006 


—-FTNsLsEe 
{ PROGRAM REVE1G2 
re DIMENSION UP(25925)3G6(25)s AUKX(255925) sUPP(25525)sFK( 25525) i, 


DIMENSION AUX1 (25925) sAUX2 (25,525) % 
CeFINO THE K MATRIX DERIVED FROM UPP=G*#UP WHERE GsB, TO THE MINUS ONE-= “ioe 
C AND K=(I+UPP)INVERSE TIMES (I-UPP) 

1 RFAD 8995N 
IF(N) 29253 
2 STOP 
3 PRINT 8989N 
898 FORMAT(16H1UP MATRIX, N=13/) 
899 FORMAT(I5) 
DO 48 I=l»N 
READ 897s(UP(I5J) sJ=l9N) 
4&8 PRINT 897s(UP(IsJ) sJ=15N) 
897 FORMAT(4E20.10) 
PRINT 8965N 
896 FORMAT(15HOG MATRIX, N=13/) 
READ 89559(G( J )sJ=lo9N) : 
47 PRINT 8959(G( J}sJ=1leN) : a -% 
895 FORMAT (4E20.10) 


C 
DO 46 I = l»sN | 
DO 46 J = 1sN 
MemUPP{1,J) = UP(I>+J) * G(T} 
PRINT 894 


DO 45 I = 1s5N a 
45 PRINT 893515 (UPP(IsJ) sJ=lsN) | 
894 FORMAT (11HOUPP MATRIX/) 
893 FORMAT (15/(1X7E17¢10)) 


DO 44 I = 1»N 
DO 44 J = 15N 
IF(I-J) 43542543 

42 AUX(IsJ)= UPP(IsJ) + le 

» AUX2(15J)=1.e-UPP(I5J) 
* GO TO 44 

43 AUX(I»J) = UPP(IsJ) 
mex2{(}.J)=UPP(IsJ) 

44 CONTINUE ; 
CALL MATALG (AUXsAUX19NoNo1l9DET925) 
CALL MAMUL (NoNoN9AUX1 9AUX29FK 925925925) 





PRINT 892 
892 FORMAT(Q9HOK MATRIX/) 
DO 40 1 = 1s5N 
4QO PRINT 89ls Is (FK(IsJ)9 J=1lsN) 
891 FORMAT (I15/(1X7E17¢«10)) 


DO 39 I=l»N 

~_ 39 G(1)=1./G(1)**2 
PRINT 8905 (G(I)y I=15N) 

890 FORMAT ( 9HOB MATRIX /(1X7F17.10)) 
GO TO 1 | 
END 
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SUBROUTINE MATALG(AsXsNR»NVsIDOsDET sNACT) 0000: 


DIMENSION A(NACT sNACT) »X(NACT sNACT) 0000 
me1DO) 1ls2ol] Y ~ 0000. 
wo 3 [=1+NR 0000 
DO 4 J=l1>+NR 0000: 
fi) 5 J)EN.0 0000! 
meri sl1)=1.0 r 0000 
NV=NR -- 900007 
hei=1.0 0000! 
NR1=NR=2 ‘i 9000: 
be S K=1>5NR1 0001! 
IR1=K+1 | oleren! 
PrVOT=0.0 : 0001. 
DO 6 I=K>sNR | 0001 
Z=ABSF(A(I9K)) 0001: 
IF (Z-PIVOT) 696597 7 0001! 
Ve = Z o0olt 
IPR=] 0001° 
CONTINUE | 0001! 
TF (PIVOT) 85998 7 | 0001' 
DET=0.0 . 0002( 
RETURN : 0002° 
IF(IPR=K) 10511510 O007, 
DO 12 J=K>5NR 0002: 
Z=A(IPRsJ) 0002: 
AC(IPRsJ)=A(KsJ) 0002! 
Park sJ)=Z | 0002¢ 
Peels J=l1.NV te 00027 
Z=X(IPRosJ) 00028 
X(ITPRsJ)=X( Kod) . 0002: 
X(KoJ)=2 ; O00 31 
DET=-DET 0003: 
DET=DET*A(K 5K) 0003: 
PTVOT=1-0/A(K5K) 0003° 
DO 14 J=IR15NR 000 3¢ 
eke) =A(K,J)*PIVOT 0003! 
DO 14 I=IR1>5NR 000 3¢ 
PHOT sJV=AlCITsJ)—AlIsK)*A(K sd) O00 3) 
DO 5 J=lsNV | 00038 
maeextKsJ)) 1595515 000 3° 
pepe) =X(KoJ)*PIVOT ; ne 0004( 
DO. 16 I=IR1>5NR . 0004° 
KT 9S) EXC 19 J)=A(T 9K) ¥X (Ko oF 0004: 
CONTINUE 0004: 
TF(A(NRsNR)) 1799917 0004 
DET=DET*A (NR 9»NR) ae | 0004: 
PIVOT=1e0/A(NR>NR) O004E 
DO 18 J=l»NV , 0004" 
X(NRoJ) =X(NRoJ) ¥PIVOT | : 0004: 
DO 18 K=l1>5NR1 0004: 
I=NR-K 0005/ 
SUM=0.0 . 7 \ 0005 
DO 19 L=I.,NR1 0005: 
SUM=SUM+A (I 5L+1)*#X(L41 9J) 0005: 
X(I9J)=X(15J)-SUM 00054 


eye 42 00055 


} 


SUBROUTINE MAMUL(MoNoL 9A9BsCoMDyNDsLD) 
DIMENSION A(MDyND) 98(NDsLD) 9C(MD5LD) 


1 S=StA( 1 sK)¥*#B(KsJ) 


2 ClIsJ)=S 
RETURN 
END 
END 
FINIS 
EXECUTE. 


4 

924067179924E-02 
321644217545E-01 
6 299614/215E-01 
726312984866E-01 
2eNNMDONNONNOE OO 


—~223.1546/15249E-01 
—6.21823595724e—07 
=—4.6 307 54722 VOE-—On 
6e03Z2Z0762750E=01 
4eQN0N0N0NTNNNE AN 


43 


3.7708598062E-01 
5 66007619310E-01 
-7.0046990430E-01 
2 .3122009925E-01 
6.NNNNNNNNANE oN 


—-8 .9024101645EF =! 
4e3819719533E= 

—-1.e2305981134¢8— 
le (61912215 0E 
BeNANNNNNNNNNE 


APPENDIX 8 
NUMERICAL RESULTS 
The following 3 x 3 unsymmetric matrix was tested in REVEIG so that 
the results could be compared with the known eigenvalues and eigenvectors. 


This example was taken from Lanczos [2]. 





The computer results compare exactly with the above answers as 
shown in the accompanying data sheets. 
The following 4 x 4 symmetric matrix was tested in REVEIG so that 


the results could be compared with its known eigenvalues and eigenvectors. 


2.8 -0.8 0.0 0.0 











An additional 3 x 3 symmetric A was tested in order to verify the 
program and the theory for 3 x 3. These results are also shown, and they 


contain the same UP and VP as in REVEIG1 for the same case. 


45 


= ae aoe oe ee 








52553 91010! 30 52567 O101022C 
627513 Q3C1004u0 52513 Q2CI1NTHO 
52513 O301N<40 57513 QO2Q0710uULUC 
52743 O7CLODRC Bo Tas GROROL Re 
62771 Q7CLOULC 52056 OTQLOCID 
52736 O7CFLOAT B27 S40 GOED IC 2 
52756 Q7TGLOLC2 7 Too memLOCe 4 
53010 Q70LDI D2 52017 Q7OLORD2 
W757 SCALE 47557 CXSOQRT 
u765G QSQXMCOF. u7N56 XTOI 
W7K55 QSCXTOL ~ 47377 LOGF 
47364 QZQLOADA 47347 OVELU 
TRACE = 6.000C0GQCE CA 
ALMOST TRIANGULAR FCRM 
TRACE = .6.00000C0GE CC 
LAGUERRE ITERATICNS 
REAL PART IMAGs PART 
MEERATE +.154U7005HE 00 a 
ITERATE 3.G0016635E dO 0 
EIGENVALUE 3,0000CCCIE 00 0 
MTERATE Dee 1 ae nae i 
[ITERATE 1.9999S99SE 00 C 
EIGENVALUE 1.999.99S99E 90 0 
ITERATE 1. COOOCOOCE CO 0 
EIGENVALUE 1.,CCOOQOQD00E 00 0 
SUM OF EIGENVALUCS = 6000000. Cgc 
@U MATRIX 


L | - 
~3.7500000002E 06-3.999999S94é GC-U.CCODOCOIOUE CO S OVTHOGONAL 
3.Q0000000007E 00 3.249999995CE 00 3.CCOOOCOIT6E 00 





1.000000000CE 00 1.000000000CE 06 1.C¢0000C0C00E 00 
| TRACE = 6.00C000C0UE CC. 


ALMOST TRIANGULAR FCRM 


TRACE = 6.0000CCOCE C9 
EIGENVALUE 3 .COCCCOOCE 00 "0 
EIGENVALUE 1.00000000E 00 0 
EIGENVALUE -- -2.00060000E 00 0 

a : - SUM OF EIGENVALUES = 6.00GCCOGCE 00 
GV MATRIX 7 


| 
2.4999999999E-01-1.N551915228E-10 1422353333216 00 

0-3.3333332352£-01 1.3333322218E 00 | — CRTHOGONAL 
ODDDCE 00 1.00C000000CE 00 1.C0000CCCO0E 90 


WO 
— J) 
«x © 


499999913E-02 8.4401108322E-10-2. 61934474 11E-09 
03850456E-10-8 2333333174 2E-02-5.E7897 27523E-09 << ORTHOGONAL 


] 
See CET 1O~ 20 O90 0 eae 





- sR ae 


SN 





nr: 
i) 


p> -COCCOSOCOCE OG-1.G000GCCOUCr GO= luc G00CCUCOE aa 
GW MATRIX fi 


C 
RI 
6.2499999S12E-02-S. UOT TOF 32 2E=-10 2.61934UN7411E-C9 
-2.9103830K56E-10 6§.3333331742F-02 5.8789737523E-C9 
4. 656612873CE-10 2.0954 75792SE-09 3. 23333325 2uE-01 
U-PRIME MATRIX : | 


-1.3000000012E O1 1.2856406521E 01 6.9282033522E 00 7 
ei. -Z000000011E O1-1.1258330335E C1-5. 1961525318 coe ORTHONORMAL 
4. OOCOO0CO28SE 00-3.4EN1016U8ZE 00-1.73205C8338E 006 

V-PRIME MATRIX 


TS. 


1.0000000007E 00-5.0N0931352EE-10 2.309%011096E 00 : 

7 6-1. 15470C5501E 00 2.2094011C91E 00 <- ORTONORMAL ! 

Le ii ooc e326: OO 2.464101648ZE 00 1.7320508238E 00 7 | 

V-PRIME-TRANSPOSE TIMES U-PRIME 


9999999995 1£-01-1.2107193N7CE-O08 147927959561E-0E : 
-~3.7252902995E-09 9.9999999951£-01 3.50u1011e70E-08 = 3 
3.2596296112E-09 1.2522100722£-08 1.CCOOOCO008E 06 - | 
K MATRIX | | 


9.26985 1986GE-02 1.198425295S£ 00 &.6279292227E-02 : 

1.0378667465E OC 5.1C645085592E-9 1-8. 98819955 36E-01 | 

| 4. 634926299 5E-02-5.992126H054UL-0} 9.5686024195E-61 
KT MATRIX : 


| 
meme O9S525495E-02 1.0278667592E CO 4leoZu9Z2oc29g9fE-02 
1. 19842529K6E O00 5.104508 746SE-01-5.9921 2647 1UE-Ol 
m2 19277416E-02=8.9681896684E-01 9.5986026iZ2Z/E-01 


SORT TG ye ee ree aha aati 


ore ree: 


TRACE = 5.69000G00E CC | 

ALMOST TRIANGULAR FCRM | 

TRACE = 5.600CCUCCE 00 | 

LAGUERRE ITERATICNS | 

REAL PART IMAG. PART ; 

[TERATE 3.4199909SE 00 G 3 

ITERATE 3519399497E 00 C ! 

EIGENVALUE 3.1937784SE 00 Co 

ITERATE POSSE Dae ae Oe C ; 

ITERATE 15616602 scene 6 | 

: ITERATE 1.611 7804CE 00 @ | 

EIGENVALUE 1.611 7804 CE. 00 6 | 

ITERATE 1.03874635E 90 0 

— ITERATE 6.85642675E-01 C : 

k 

EIGENVALUE 6 .85642675E-01 0 ; 
ITERATE | | 1.08798451E-01 ; 0 


EIGENVALUE — 1.879845 1E-01 “6 








SU! UR se 1uUENV AAszwEe SD = PeOVUUUUGUYE wu | 


8U MATRIX 


OG=5..052697 I049 nec: 
00 2.48705464H2TC Cl 


1 .2326497262E-01-3.92F0S8177CE=01 146303529485E 
4. 146636067 2E-01-1.040807769ZE 00 2.4222642054HE 
7 «2800 385559£-01-7.1410660635C-01-3.C294SCTSSIE 00-6. 984UN62315C 00 
1. 000000000CE 00 1.00G000000CE 00 1.CCOO0CCOOOE 00 1.C000000000E CC | 
TRACE = 5.600C00G0E 09 | 

ALMOST TRIANGULAR FORM 
5.600000fCE CO 








TRACE = 
| LAGUERRE ITERATICNS 

REAL PART IMAG. PART : 
ITERATE 2.419900S9E 90 | 
ITERATE 3.19399N97E 00 0 : 
EIGENVALUE 3.193778URE 06 C | 
ITERATE 2.45552682E 00 @ | 
[ITERATE 1.61669232E 00 0 ! 
ITERATE 1.613 78940E 00 0 : 
CIGENVALLE 1.61178040E 00 0 3 
: 
ITERATE 1.03874635E 00 G ! 
ITERATE 6 .85642675E-G1 0 | 
Bee IGENVALUE 6.85642675E-01 
ITERATE 120879845 1E-01 0 : 
: EIGENVALUE 1 .08799%51E-01 C 
| i 
SUM OF EIGENVALUES = 5.6Q00COUCCE CO | 
BV MATRIX : 
] | 

) 1.2326497265E-01-3.93E058i1770£-01 1.62C852948SE CC-5.C526979089E Ci 

He 14663 00672E-O1- 1.040807 7692E CO 2.4222643G5KE 00 2.487054N6R2TE C1 

7 .2800385559E-C1-7.141066663£5E-01-3.029U507S51£ OC-&.984N4E2515= CC 

1.000 9000000%E" 00 1,.000000000CE 00 1.00000CCCOUE 00 1.00C0CCOCCOE 06 





4 


Mer if i297 7KOE OO 2.67H6420I19CE-0O8 6. f2462656/77E-O8-1.0K 77578964 E=C>9 
26746420 190E-08 2.7483121654E 00 9.26665961 THE-08-7. 5669959188E-69 
Reg 02 05677 E-08 Fe 2OCOS9OITHE-CE WeS7OR6IISS82E 01-6.519258022sE2e3 
meme? (317 8964E-09-7.5669959 1S8E-99-6.5192580225E-C8 2.221502! Gabe Oe 


U-PRIME MATRIX i to 









3.7708 598062E-01-&. 9024101 6uSE-CI 


1 
9.406717 992ZUE-02-2.375467 S24SL-01 | 
366076 19310E- Odie 3819719533E- -C1 


3. 164421754 $E-01-6. 278233 9352 4E-01 
(565556 147215E-01-4. 207547 221 5t-01-7, 0084990430E-01-1.2305981130E-C1 
7.831 290N866E-01 66052078294CE-01 262122009925E-01 | 167619122150E-02 | 


J-PRIME MATRIX 2 
3 | 


| 9.806717992NE-02-2.37S4E7524SE-O1 3.77C08598C62E-01-&. 9024 101645E-01 
| 3. 1648217545E-01-6.2782339524E-01 5.6007619310E-01 4.38197 19523E-01 


fa ras. = S50) = a > 20 =—({)j— 2 0) & 3 OE-C ] 


« 










Wer RIME-TRANSPOSeE TIMES U-PRI-E 


1.QOO000000CE 00 1.23N002H1TNE-G8 1.18839 14742E-CB-1.0913536U21E-11 
1.2340024 11LE-08 9.95995999955-O} 1.29330 1NE59E-08-7.95807864C6E+ 11 
1. 1888914 742E-08 1.293301465SE-08 1.COOD00GCCOE 00-2.6886937121E-I1C 
-1.0913936421E-11-7.958078640 CE-1 1-2. 6886927121E-1C 1.COoO0000OCOE CC 

-. -K MATRIX Oe | | q 


~.9017760876E-09-1.499670178EF 00 3.3570411263E-01 1.56168075S1E CO 

I.¥9967TOIZIZE 00 1.6872831489SE~G8-4.G592364591E 0G 3.9046217502E-C1 
“3. S97ONTIUUNE-O1 4605923652126 00 3.2716343N12E-CE-1.9242020291E CO 
~1.5616867415E 00-3.9046223332E-01 1.9242020085E 0G 1.86992110E8E-CE 
KT MATRIX | 


4.0017766872E-09 14499670173€E 90-3«5970411263E-01-1.5616807591L CC 
-~1.8996701212E 00-16872831NB9SE-08 4.05923EN591E CO-2.90462176025-C1 
3.397041 144 E~01-4.0592365214E 00-3.2716343412E-08 1.9242020291& CC 
1.5616807415E 00 3.9046223332E-01-1.9242020085E CO-1.8699211068E-Cea 











OVUM UF CIVECNVAaALYUEDS =- Z~eevuyuvuvuuvuuvunc wwe 


U MATRIX 


1 
3e-O0772N17209E-01-5.822811950CE-01 7.4¥I2236765E 00 
ONLI a ee eC) 1 ilo MIE oI 00 
1T~-OQ00Q0000000E 00 1.000000000CE 00 1.0000000000E O00 


TRACE 


= 5.2000000E 00 


.. | ALMOST TRIANGULAR RORM 
TRACE = 5.20000000E 00 
LAGUERRE_ITERATIGNS 


REAL PART IMAG.. PART 
[ITERATE 3.32218337E 00 0 
ITERATE 3.19296062E 00 0 
EIGENVALUE 3.19293436E 00 0 
[TERATE 2.16568534E 00 ) 
ITERATE 1.547 13834%E 00 ‘oe 
EIGENVALUE 1.54713834E 00 
ITERATE 7 +69321492E-01 0 
ITERATE 4 .59927299E-01 0 
EIGENVALUE 4 .59927299E-01 0 

SUM OF EIGENVALUES = 5.20000000E 00 
V MATRIX 


3.07724 17209E-01-5 8228 119506E-01 7e4H12236765E 00 
9.0012116814E-01-9. 118972315 ZE-0 1-3. €548905973E 00 
1.8000000000E 00 1.0000000000E 00 1.0000000000E 00 
W MATRIX 


1.9049122834E 00 2.91038304S7E-11 133877620 10E-09 
2.9103830457E-11 2.170607951CE 00 2.0954757929E-09 
1¢3387762010E-09 2.0954757929E-09 6.9730035083E 01 
U-PRIME MATRIX i. | 


s 


Pad 


2. 2295854460E-01-3.952228721K4E-01 8.9111636088E-01 
6.5217400856E-0 1-6. 1894 94 8008E-0 1-4. 3768779841E-01 
%.2u5u023707E-01 6.787491601SE-01 1419754008 16E-01 

V-PRIME MATRIX | 


2 .2295854460E-01-3.9522287214E-01 8.9111636088E-01 
6.52174 004S6E-01-6. 1894948008E-0 1-4. 37687798U1E-01 
7. 2454023T707E-O1 64787491601SE-01 1.19754,00816E-01 
V-PRIME-TRANSPOSE .TIMES U-PRIME 


1 
9.999999999E6E-01 3.637978807 1E-11 1.1641532182E-10 
3637978807 1E-11 T.OOOOHODOOCE OO 1.74622982 THE-10 


1.1641532182E-10 1-746229827HE-10 1.00000CCOOOE 00 
ee 


a) 


YP=UP 


= I. 








ee eee cai saan! 


. 





(-1.27329258256-11 LY 7ISR3H72ZE 00=-2.3015283608E-0) 
~TeURTISHZUTZE 00 2.764863893272E-10 1.542544 7924E 00 
2. 3015283620E-01-1.54254N7925E 00-1.6734702512E-10 

KT MATRIX 





1 
1262732925825 E-11-1.44715U3U72E 00 2.3015283608E-01 
Penge VON on Ze 00-2.7648 63893 2E-10-1. SU254N792KE O00 
—2.3015283620E-01 1.54254U4U7925E 00 1.6734702512E+10 


te Oy 


The K matrix derived from REVEIG using a 4 x 4 symmetric A was used 
as input into REVEIG1. The results for UP and VP from REVEIG] compared 
with those of. REVEIG. In addition, UP equalled VP as required in this 
case. 


Other sample results are included, as discussed in Section 4 of the 


text Qa 


a2 











K MATRIX N= h oe SKEW? SKM ET 
~K.OO01TTO687TBE-09 == 1.49967.-'178 
1.499670121ZE 00 1.872821489 
-3.397041] JUNLE-01 4.0592 4652 1: 
-1.561 68074 15E 00 ~3.904622338 
U-PRIME MATRIX = V PRIME 
1 


9.4067179903E-02-2.375467525Z2E-0 1 


3.1644217541E-01-6.2782339526E-01 
5655561472106 01-4. 4075473212E-01- 
7.6312984867E-01 6.0320782942E-01 
V-PRIME MATRIX = U PRINE 


9.4067178376E-02-2.3754675862E-01 
3.164421 7654 E-01-6.2782340635E-01 


EC 

E vO 3.39/0411263E-C) le oot 
E-08  -4.0592364591E 00 3.9C4 
E-O] 1.9242020085E 00 1.86S 


3.77085 98059E~01-8.9024101 645E-Q) 
5.6007619314E-01 4.3519719532E-G1 
-7.00469904 33E-01-1.2305981130E-01 
2.3122009924E-01 1.7619122184E-02 


3.7708598232E-01-8.9024101 633E-C1 
5.6007619757E-01 4.3819719545E-01 


DS NR a a SAU LCS Sal) Saal CIO Lol Cet CAS Ss Jorea EIEN GIs 18 


fees 1 296S5852E-0!1 6.052078 1 7O0SE-0} 
VPT X VP 


2235122008239E-01 1.7619122242E-02 


| | 
1-Q00000000iE 00-1.2267264538E-08-1.1859810911E-08 8.1854523157E-12 
~1.2267264535E-08 1.000000000CE 00-1.2958480511E-C8 3.81987 774 74KE-) I 


-1.1859810911E-08-1.29584805i1E-08 
8.185452315S7E“12 3.819877 7474-11 
VP X VPT 


} 
1,.0000000018E 00 4.329194 7804E-09 
4. 329194 780NE-09 1.000000009SE 00 


5 2 3296389524E-10 2.2491803975E-09 


=4.4537955546E-09-4 «B44 878276 2E-09 
UPT X UP 


OOOOOOOOTE OO }. 236912794NE-08 
369127944 E-08 1.000000000CE 00 
97965271E-08 1.2914824765E-08 


— oS 


oa 
-~- gy YW —-— — 


99999981 1E- 01-4. 336470738 1E-09- 
364707 387E-09 9.999999901 3E-01- 


e 0 - 
WWNHDO x OF ~WNNO- 


-5,8933u79986E- 10-2.3073880584E-09 9.9999999269E-0 
424110493036E-09 4.8153196985E-09- 
VPT X UP 


1 


1.0000000007E 00 4.3655745685E-11- 


5 -8207660913E-11 
~4.7293724492 E~11- 30274 1809264E-11 


17155614E-12-4.092726 |S8CE~1 I- 


9.9999999993E-01 2.5011104298E-10 
2.5011 104296E-10 9.9999999998E-01 


2296 389523E-10-4.4537955546E-C9 
2.2491803975E-09-4. 8448782 763E-C9 
1.0000000073E 00 1.1539555089E-08 
1.1539555089E-08 $.9999998095E-C1 


NX 


1.1797965271E-08 5.911715561KE=12 
1229148247 65E-08-4 .0927261580E-11 
1.0000000001E 00-2.1941559680E-10 
2.1941559680E-10 1.0000000000E 00 


-5.4933479986E-10 4.4116493036E-09 
2. 3073880584E-09 4.8153196985E-09 
1-1.1431325220E-08 
1.1431325220E-08 1.0000000192E 00 


1.0913936421E-11 9-O949L7OI75E-12 


1T-0O00000000CE 00-2.5465851650E-11 3.774H03012KE-11) 


1.0CO00000000E 00 5.6843418860E-12 


8 .5401996666E-12-3.092281986CeE-1 I 2.5011 1O4299E-11 9.9999999998E~0 } 


53. . a! 
Py — ont * — oe ty at + omen t wre rey wee oe nt 





2.2295854460E-01-3.9522287215E-01 8.9111636085E-01 
6+5217400455E-01-6.189494801 3E-0 1-4. 3768779843E-01 
7.2454023702E-Q1 6.7874916019E-01 1.1975400812E-01 
V-PRIME MATRIX & U PRIME : | 


24229585445 1€-01-3.9522287 22 SE-01 8.9111636090E-01 
66521 7HOO0463E-01-6. 189494 7997E-0 1-4. 37687798 39E-01 
Te2454023706E-01 6.7874916012E-01 1.1975400796E-01 
VPT X VP 


1.0000000000E 00-3.455191522EE~-1 1-1. 27 32925825E-10 
—-1.4551915228E-11 9.9999999988E-0 1-1.4006218407E-10 
~1.3732925825E-10-1464006218407E-10 1.00000CO000E 00 
VP X VPT . 


1'.0000000000E 00-2.182787284 ZE-1 1-1. 2369 127944E-10 
-2.1827872845E-11 9.9999999998E-01 1.3733369996E-10 
-1.236912794K4E-10 1.3733369996E-10 9.9999999991E-01 
| UPT X UP ~ 


9.9999999991E-01-249 10383045 7E-11 Te8216544S552E-1)1 
-2.9103830457E-11 1.-0000000001£ 00 1.7462298274E-10 
TeB21654N552E-11 1.7462298274E-10 9.9999999993E-01 
UP X UPT 


9..9999999988E-01 72275957614 1E-12 1.0913936421E-10 | 
7.2759576141E-12 1.0000000001E 00-1.6916601453E-10 
1.0913936421E-10-1.6916601452E-10 1.00000CGCO00E 00 
VPT X UP ~~ 


: 
9-9999999995E-01-3. 637978867 1E-11-3.4560798667E-11 
(776275957 61UIE-12 9.999999999FE-01 5.4569 682105E-12 
7. 275957TOIGIE-12 2.9103830457E-11 9.9999999998E-01 


1 
1H 


‘ 











—— a oo 














SKE WHSTMMETIUIC 


K MATRIX, N= = 
-  Q -1.0000COCCOOE 00 -4.0000000CO00E GO -5.0000CO0CC 
1. COOOCOOL00E 00 GC -3.Q0000000COCE 00 -4.CCCOC0000 
4 OCOOGOOLOCE 00 3. OOOICOCOOODE 00 0 -3.000000000 
5.000GOQ00000E 00 sOOOOCCOOOOE CO 3.COO0000000E 00 
i 
U-PRIME MATRIX Orie 


=e 99002469 15E-0 1-6. 1728395055E-0 1-2. 9506 172836E-01 
—9.6296296294E-01 2.5925925934L-0 1-7 -UCTHOTHIO2ZE-02 
Jes mee OOF} rE -01-7.1604938272E-01 6.1728394998E-02-6 .6666666665E-C1 

l. 2345679009E- -O1 1.975308642Z2E-0 1-9. 1358024699E-0 I- 


V-PRIME MATRIX B® U PRULME 


626666666665E-0} 
0 


2,33355555500E= 00 


~ 1. 35802469 15E-01-6«172839505£E-01-3. 9506172836E-01 6 -6666666665E-01 
—9.6296296294KE-01 2.5925925934E-01-7-4O7HOTYHIOZE-O02 
-149753086417E-01-7.1604938272E-O1 °6.1728394998E-02- 6. 6666666665E-01 

T.23545679009E-0} 16975308642ZE-01-9. 1358024699E-01-3.3333333330E-01 


VPT X VP 


0 


1 
9.99999999935E-01-9. 458744898 4E-11 5.4569682106E-11-1.7280399334E-11 
—9.4U5874UU89BY4E-11 9.999999999 BE-0 1-4. 729372UN92ZE-11 4.1856756282E-1 1 


524569682 106E-11-4. 729372449 ZE-1 1 


1-000000000TE O00 


6.54835618528E-1 1 


-1.7280399334E-11 4.1836756282E-11 6.5483618528E-11 9.9999999995E-01 


VP X VPT 


1 - = 
9.9999999988E-01-1.6 37090463 2E—1 1-7. 275957614 IE-12 3.6379788071E-11 
-1.6370908632E-11 1-000000000CE 00-8.6060926155E-11 &.5492501967E-11 
—7.2759576141E-12-86 6060936155E-11 9.9999999998E-017 1.4551915228E-11 


4 
3.-6379788071E-11 8.5492501967E-1 1 


UPT X UP 


1264551915228E-11 


1.-Q00000000I1E CC 


9. 9999999993E-01-9.4587448984E-11 5 64569682106E-11-1.7280399334E-11 
-9 458744898 E-11 9.9999999998E-0 1-4. 72937 24U92E-11 4.1836756282E-11 
3° 4569682 1O6E- W- 4. 72937T24UU9ZE-11 1.-000000000TE 00 6.5483618528E-11 
i. 7280399334 E- 1 4.183675628ZE-11 6.5483618528E-11 9$.9999999995E-0 1) 


UP X UPT 


1 ° ‘ : 
9.9999999988E-01-1.6370904632E-11-7.27595761KIE-12 3.6379788C071E-1 1 
ae umes Oe E— I 1-0000000000E 00-8.6060936155E-11 8-5492501967E-11 
ag fae ON TE 12-8. 6060936155E-11 9. 99999 979 38E- OF 124551915228E-11 


3. ear CCOTIE~ 11 8. 5492501967E- 11 


VPT X UP 


le 45519,15228E- 11 


1.0000000001E OC 


9. 9999999992E-01-9.45874489BNE-11 5.4569682106E- 11-1. 7280399334E-11 
-9.4587448984E-11 9.9999999998E-0 1-4. 72937 24492E-11 W. 1836756282E-11] — 
5.569662 106E-11—-4 6 729372449 2E-11 1.0000000001E 00 6.5483618528E-11 
~1.728039933UE-11 4.1836756282E-11 6.54836 18528E-1]1 $.9999999995E-01 


55 ee 





MATRIX, N= 4 mUNSYMMETEIC | 
9.0000000000E O00 2-0000000000E O00 5.0000000000E 00 1.000000000: 
1.0000000000E 00 6e-0000COO0000E O00 2.Q00Q0000000E O00 2.0000C0000; 
&.OOD00000000E O00 3.0000COO0000E O00 TeOOQ0000000E 00 3.000000000. 
5. 0000000000E OO 4 ,OOD0D0000000E GO 3.0000000CCOE 00 4 .0C0000000; 


U-PRIME MATRIX 


5017921147 5E-02-6. 315412 186 3E-0 1-8. 3154 121875E-02-1.07526881 72E 701 2 
3 


-4. 3010752691 E-02-3.0107526872E-02-6. 43010752 65€- =(1e o354030712E-$1 
u 
~2, 50896057 3WE-01-2.4229390682E-01 1. 5770609340E-02—4 «6236559 140E-01 


V-PRIME MATRIX 


1 | 
~7.6344086018E-01-3.4408602157E-02-1. 63440860 19E+01 644516129030E-02 
2 i 
| | | : 
—-1.2068965517E 00-2.0689655 17 SE-0 1-1-3793 103452E-01 7.5862068973E-01 
2 ‘ 
1. 7733990 1U9E-01-126798029557E 00-1. 67487684 72E-01 7.7832512329E-01 
3 : 
2.75 8620689UE-01 2.7586206902E-01-1-4827586206E 00-3-4482758626E-0 1 


: -3.2512315277E-01 2.463054 1874E-01 6.4039408874E-01-2.0935960592E 00 
_ vPT x vP i 

| : 

1.6698536727E 00-5.2173068997E-02-4. 8047 7565S4E-01-1.9199689376E-01 | 

: ~5,2173068997€-02 3.0013103936E 00 5.8579436492E-02-2.0751777529E 00 : 


3 

“4, BON7TSOSSNE-O1 5.8579436U92E-02 2. 65575¥B108E 00-1. 0644276738E 00 
21..9199689376E-01-2.0751777529F 00-1.064%276738E 00 5. 6833458720 OC : 
VP X VPT 


1 : | 
2.0939357908E 00 76470698 1507E-0 1-4 647086801 4HE-01-1-3351452355E 00 : 


2 

7.4706981507E-01 3.4870295327E 00-4. 3451673200E-01-2.2081584125E 00 
3 . a 

~4'.4708680 1ULE-01—4e 345167 320CE-01 2.4696789536E 00-2.4936300312E-Cl 


uy 
—1.33514852355E 00-262081584125E 00-2 .4936300312E-01 4.9596204719E 00 


| Me 15865669E-01 S.6464225799E-02 1.4430441538E-01 6.9680502568E-02 
; Pe 5225799E-02 4 .596410631SE-91 7+3677367961E-02 1.8354337690E-01 
a Me s0uyts3sc-01 7. 367736796 1E-0Z 4.4733906290E-01 1.1555863876E-01 

6.96805025686-02 1.8354337690E-01 16 1555863876E-01 2.6696727946E-01 


| 
| , 
UP X UPT | : | : 


| 
| UPT xX UP 





ee a 


F 6-61490114459E-01-9.9248468005C-03 1.2647936175E-01 1.6747369637E-01 
7 2 
-9.9248468005E-03 4.1983890238E-01 9.1136547578E-02 1.8883416196E-01 


3 
1.2647936175E-01 9.113654U757BE-02 453680194 20E-01 9.7435541679E-02 
4 
1.6747369637E-0O1 1.8883416196E-O1 9.743554 1679E-02 3.3568582111E-01 
VPT X UP 


1 
9.9999999995E-O1 2.3646862246E-11-3.0468072509E-11 2.91058304K57E~-1 1 


2 : | 
=] BIBIBINOSSE- 12 }-000000000CE 00-3. 29691 82939E-12-2. 728484 1053E-11 = 
3 = 


1. isieumns Gl 11-5. 4569682 106E- 11 9.9999999993E-01-2. 9103830457541 1 


a): 4551915220€- Tiee 910383045 7E-T I~ ~3.6379788C71E-11 1.-0000000001E 00 


t 








N= gy SUNSYMMETIC WITst ZERO PLAGONAL 


K MATRIX 
‘ 0 2. 0000000000E CO 5. 00000000C0E 00 1.0C0000000 
1.0000000000E 00 0 2. OOO0000CO0E 00 2. CC6CC00000 
4 OOOGOOOCOOE 00 3.0000C00000E 00 0 3° 000000000 
“O000000000E 00 h -OOOOCOO000E 00 3.0000000000E 00 
U-PRIME MATRIX | t 


f 


-3.4999999998E 00 4.2499999997E 00-2.7499999998E 00 2.2499999999E 0 

3.4999999998E 00-7.750000000CE 00 &.2499999999E 00-2.7499999999E 00 
-4.9999999964E-01 127499999995E 00-2.2499999998E 00 7.4999999986E-41 
~1.291N824766E-10 5.000000001CE-01 4.9999999990E-01-1.4999999997E 00 
V-PRIME MATRIX | 


~6  4666666666E-01-3. 333333333 TE-O1- 1. 6666666667E-01~1.6666666668E-D1 
~3. 7500000001 E-01-3~7500000002E-0 1-2. 8124999999E-0 1-2. 1875000001E-01 
4. 166666667 1E-02-2.9 16666666 9E-0 1-8.0208333333E-01-3 .6458333334E401 
=229166666567E-01 4. 1666666656E-02-143541666668E-01-6.9791666666E-01 
vPT xX VP 


6.7187500002E-01 3+3854166670E-01 242265625001E-01 3.8151041669E-01 
3, 3854166670E-01 36385%166673E-01 3.8932291670E-01 2.1484375004E-01 
2+2265625001E-01 3.893229 1670E-01 7.6855468751E-01 4.7623697920E-01 
3.815104 1669E-O1 2¢1484375004NE-01 &.7623697920E-01 6.9563802086E-01 
VP X VPT 


6. 1V11111115E-01 4. 5833333338E-01 2.6388888892E-01 3. 194RRkuky 7E-CI 

4.58 33333338E-01 4.0820312502E-01 3.9908854168E-01 2.8450520835E-01 

2-6388888892E-0 | 3.9908854168E-01 8.6306423615E-01 3.3875868057E-C} 

3. T9UULUEEYETE-O1 2.8450520835E-01 3.3875868057E-01 5.9223090279E-01 
UPT X UP 


2 .U7%9999996E 01-4. 2874999995E 01 265624999997E 01-1.7874999997E 01 
~4.5874999995E 01 8.1437%99992E 01-4.8312499996E 01 3.1437499997E 01 
2 .3624999997E 01-4.8312499996E 01 3.0937499998E 01-2.0312499997E 01 
~1.7874999997E 01 3.1437499997E 01-2.0312499997E 01 1.5437499997E 01 

UP, X UPT 3 oe re : 


| We2937K9999NE 01-6.3062K9999NE 01 1.706249999KE 01-2.6249999979E 00 
 =6. 306249999KE 01 957937499997E 01-2. 6937499993E 01 2.3749999970E 00 
1.7062499994E 01-2.6937499992E 01 8.9374999970E 00-1.3749999992E CO 
-2.6249999979E 00 2.37N999997CE 00-1.3749999992E 00 2.7499999990E 00 
VPT X UP . 


99999999991 E-01 2 58296N'953CE-10-5.8207660913E-11-2.9103830457E-11 
-7.550277%320E-11 140000000002E 00-1.9099388737E-11 2.3646862246E-11 
~2.2261777606E-10 366379788070E-10 9.9999999993E-01 6.5483618528E-11 
-2.2642628509E-11 1.7462298274E-10 3.6379788071E-11 9.9999999979E-01 


{1 
ra 














Only one test is shown for REVEIG2 since the results present no 
particular significance except that the matrix K is not skew-symmetric 


as was anticipated. 


>) 











G MATRIX, 
2-QOCOOCOOQO00E OG 
UPP MATRIX 


Pa 
tt 
= 


1. 
Ve 
Bere > 50C8S29E 00- 
Oe 


FUWONMNO— 


K MATRIX 


4.QOOOCCODOOE CO 


6-CO0000COCOOE O00 


8.CQ00COO0O0O0C 


B8I34U55985E-01-4.7509350N97E-01 7. 5417196124E~O1- 1.7804820529E QC 
657687018E 00-2.511293581CE 00 2.2h030K7724E 00 1.7527887814KE gc 
©584U528392SE CO-uU.ZO0281S4258E 00-7. 
1050387892E 00 4.8256626357E 00 PEROT OCTINOE 00 1.4095297720E-C1 


S8S55001 8 Cee Q1 


1 | 
Oe(4osub 712 -01 2.99332496SCE-01 &. 78359027 INE-01-2.2909850124E-C2 
ue oe E-0e 1.07997842] 4E-01-7. 18050 15088E-01-2.7720428801E-C2 


ae oo Cll =O02 8.4 (606687 1 CE-O l= 192596555 5E—Ul] 


2.139776674UHE-04 


1.5098146321E-01 7.96808K759¢E-01 2.72158C9627E-01 $.9224176882E-01 


B MATRIX 


2.4999999999E-01 &.Z50000000CE-G2 2. 77777777 77E-02 1.5625000000E-G2 


60 























